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^ . ABSTRACT 

i-C ■ The Chandra X-ray Observatory has revealed X-ray bubbles in the intracluster 

medium (ICM) of many nearby cooling flow clusters. The bubbles trace feedback that 
q ■ is thought to couple the central active galactic nucleus (AGN) to the ICM, helping 

?-h " to stabilize cooling flows and govern the evolution of massive galaxies. However, the 

c/3 . prevalence and duty cycle of such AGN outbursts is not well understood. To this end, 

we study how cooling is balanced by bubble heating for complete samples of clusters 
(the Brightest 55 clusters of galaxies, hereafter B55, and the Highest X-ray FLUx 
Galaxy Cluster Sample, HIFLUGCS). We find that the radio luminosity of the central 
galaxy only exceeds 2.5 x 10 30 erg s _1 Hz -1 in cooling flow clusters. This result im- 
| plies a connection between the central radio source and the ICM, as expected if AGN 

feedback is operating. Additionally, we find a duty cycle for radio mode feedback, the 
fraction of time that a system possesses bubbles inflated by its central radio source, 
O i of > 69% for B55 and > 63% for HIFLUGCS. These duty cycles are lower limits since 

■ some bubbles are likely missed in existing images. We used simulations to constrain 

the bubble power that might be present and remain undetected in the cooling flow 

fN| . systems without detected bubbles. Among theses systems, almost all could have signif- 

icant bubble power. Therefore, our results imply that the duty cycle of AGN outbursts 
with the potential to heat the gas significantly in cooling flow clusters is at least 60 
per cent and could approach 100 per cent. 
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1 INTRODUCTION 



A much debated topic in extragalactic astronomy in the last 
20 years is the "cooling flow" problem (Fabianlll994l b Since 
the cooling time of the ICM in t he cores of many galaxy 
clusters is shorter than a few Gyrs l|Edge et al.ll 19921 ). in the 
absence of heating the gas is expe cted to cool d own below 
X-ray temperatures and condense (|FabiarJll994T ). However, 
searches for evidence of the high mass depositions rates ex- 
pecte d in the standard cooling flow m odel were unsuccessful 
(e.g. jMcNamara fc Q'Connerj|l989h . Data for such clusters 
from the Chandra and XMM-Newton X-ray observatories do 
not show the signat ures of gas cooling below ~ 2 keV (e.g., 
IPeterson et alfe oOlb Therefore, it appears that only a small 
fraction of the expected cooling is occurring in most cooling 



flows, suggesting some heating mechanism is operating in 
these systems. 

A possible source of heating has been identified: high 
resolution images from the Chandra X-ray Observatory have 
revealed X-ray bubbles in the ICM of many nearby clus- 
ters, groups and ellipticals (very few systems were previously 
known to host bubbles in their atmospheres, e.g, Perseus, 
Cygn us A). No table examples includ e Perseus (|Fabian et al.l 
l200(t>. A2052 jBlanton et al.l l200ll ) and MS 0735.6+7421 
( McNamara et al.ll2005l b among many others that have been 
discovered. In a study of a large sample of such bubble sys- 
tems, it was found that a central radio source was present 
in ev ery system, with th e radio plasma often filling the bub- 
bles (|Birzan et al.l [2008b The X-ray bubbles are therefore 
interpreted as regions where the X-ray emitting gas has 
been displaced by radio plasma produced by energetic out- 
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flows from the central AGN and show that the AGN is 
injecting significant power into the ICM. The energy that 
the bubbles inje ct into the ICM may b e important for bal- 
ancing cooling dMcNamara et ail |2000| ; iFabian et ail 120001: 



iBlanton et al 



2001 



) and inhibit ing star formation in mas- 
Croton et ai1l2006l ). 

' Birzan et al.l 



sive galaxies 

A numbe r of authors 

Dunn fc Fabianl 12004 iDunn et al 



2004; 

20051 ; iRaffertv et all 



2006) have analyzed samples of cooling flow clusters with 



visible bubbles in their envir o nmen ts. Using a sample of 
such systems, IRaffertv et al, I l|2006h found that at least 
50% of the systems have enough energy in their bubbles 
to balance cooling, considering only the enthalpy, although 
the inte grated energy is sufficient t o offset cooling in all 
systems IjMcNamara fc Nulsenll2007l ). The enthalpy is only 
a lower limit for the total energy injected, as Chandra 
images reveal evidence for the existence of weak shocks 
dMcNamara et al.l 120051 ; IFabian et all 120061 : iForman etal] 
120071 ; IWise et al.l 120071 ) with energies comparable to the 
bubble enthalpies. Additionally, the energy may still be 
underesti mated due to further adiabatic losses, undetected 
bubbles dNusser et ai1l2006l; iBinnev et al.l [2007h or cosmic 
ray losses (|Mathews fc Brighentill2008l ). 

However, it is not clear whether the bubbles are always 
present when heating is required. Currently, the bubbles are 
primarily de tected by eye. Some of them have bright rims 
(e.g., A2052; IBlanton et all 120011 ) which make them "clear" 
cases. However, the det ectability of a bubble depends on 
its location, orientation fenfilin fc Heinzll2002l ; iDiehl et aJj 
120081 ; iBriiggen et all 12009 )7 its angular size and the depth 
of the observation. A majority of observed bubbles come 
from systems where the angle betwee n the jet axis and the 
line of sight is betw een 45° and 90° l|Enfilin fc Heinzll2002l ; 
IBriiggen et ail 2009 ) . As a result, we are likely missing some 
bubbles in the existing images. 

Indeed, many cooling flow systems, which presumably 
require heating, do not show bu bbles in their existing Chan- 
dra images (e.g, A1650, A2244; iDonahue et all 12005). It is 
therefore important to understand the biases and selection 
effects in the detectability of current X-ray bubble samples. 
However, images of many of these systems are shallow. In 
these systems, the question is whether there might exist 
enough undetected bubble power to balance heating. Also, 
there might be alternative scenarios for the lack of bubbles 
in some cooling flow systems, such as the presence of ano ther 
heating mechanism (e.g., "sloshing" ; IZuHone et al.1l2010l ). or 
perhaps because some systems are in a cooling stage, with 
high rates of star for mation, and might not need heating 
right now (e.g., A1068; iMcNamara et a71l2004l ). 

These issues are important in understanding the duty 
cycle of AGN heating, which we define to be the fraction of 
time that a system possesses bubbles inflated by its central 
radio source (and hence shows clear evidence of AGN heat- 
ing). The duty cycle gives insight into the heating process 
and is useful in constraining (or as an input to) simulations 
of galaxy and cluster formation and evolution that include 
AGN heating. 

Previous work to understand the duty cycle of AGN 
heating in clusters has been done bv lDunn et all (|2006l ). who 
found an AGN duty cycle of 70% for cooling flows in the B55 
sample. However, their selection of cooling flow systems is 
based on ROSAT data, which are less sensitive to the pres- 



ence of a cooling flow than higher-resolut ion data. For ellip - 
ticals the results are more controversial. iDong et al.1 (|20ld ) 
found a duty cycle for AGN feedback of 50% for a sample of 
51 elli pticals. This duty cycle is similar to the IDunn et al.l 
(|2010h finding of 9 out of 18 ell i ptical s having bubbles, but 
much higher than lNulsen et all l|2009l) finding by using 104 
ellipticals from Chandra archive, where only 1/4 have bub- 
bles (note however that Nulsen et al. found that the duty 
cycle varies with the X-ray luminosity, so the above results 
may be consistent). In this paper we address the question of 
duty cycle of AGN feedback in complete subsamples of cool- 
ing flow clusters using simulations to place limits on AGN 
heating for the clusters which do not show bubbles in the 
existing Chandra images. For this purp ose we chose the B55 
jEdge et al.lll990l ) and HIFLUGCS jReiprich fc Bohringerj 
2002) samples, since both samples have already been ob- 
served with Chandra. We assume Ho — 70 km s _1 Mpc _1 , 
Qa = 0.7, and Qm = 0.3 throughout. 



2 SAMPLE DEFINITION AND PROPERTIES 
2.1 Complete Samples 

In this work we pres ent an analysis o f two complete sam- 
ples: the B55 sample llEdge et al.lll990l) and the HIFLUGCS 
sample (|Reiprich fc Bohringerl 20021 ). The B55 sample is a 
2-10 keV flux limited sample of 55 systems (or 56 since we 
count both A3391 and A3395s) sel ected using data from 
the HEAO-1 and Ariel- V surveys (|Edge et alill990l ). HI- 
FLUGCS is a 0.1-2.4 keV flux limited sample of 64 systems 
based on ROSAT and ASCA observations at galactic lati- 
tude b > 20° . The two complete samples have 44 systems in 
common. For the systems that belong only to the B55 sam- 
ple, the majority of them are low latitude systems (b < 20°). 
These are: AWM7, Perseus, Ophiuchus, A2319, Triungulum 
Aus, 3C129.1, PKS 0745-191, Cygnus A and A644. M87 is 
also present only in the B55 sample, since it was excluded 
from HIFLUGCS based on the fact that the irregular X-ray 
background and the extended X-ray emission at low lati- 
tudes makes the detection of clusters more difficult. A1689 
and A3532 are only in the B55 sample due to differences in 
the X-ray band used to construct the samples. The systems 
that belong only to the HIFLUGCS sample are in general 
smaller, cooler ones (e.g., NGC 507, NGC 1399, NGC 1550, 
EXO 0423.4-0840, MKW 4 etc), and softer-spectrum sys- 
tem s (e.g., Ser s ic 159 /03, A2589, A2657 etc). 

ISun et al.1 (|2007l ) showed that there is a class of sys- 
tems having truncated atmospheres of low entropy gas, with 
short cooling time, embedded in much more extensive at- 
mospheres of higher entropy gas. Some of these systems are 
present in our sample and we refer to them as corona sys- 
tems (they are marked with an asterisk in Table [T] Tabled 
and Table 0. 

The B55 sample was studied by IDunn et al.l (|2006l ) 
in order to determine the fraction of cooling flow clusters 
with bubbles. In order to separate the cooling fl ow cl usters 
from the non-cooling flow clusters, IDunn et al.l l|2006T ) used 
ROSAT data. They defined a cluster to be a cooling flow 
cluster if t coo i < 3 x 10 9 yr in the core and Toutor/Tccntor > 2. 
They found th at 14 out of 20 coolin g flow clusters have clear 
bubbles. Also, iRossetti fc Molendil l|2010l ) studied the B55 
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sample and separated cooling flow from the non-cooling flow 
clusters using the pse udo-entropy ratio . The HIFL UGCS 
sample was studied by ISanderson et al.l l|2006l . I2009D . who 
used a temp erature ratio in o rder to separate the cooling 
flow clusters; IChen et al. I (l2007h. who used a mass deposition 
rate ratio; and iHudson et aTT^OlOh . who used the cooling 
time. However, different definitions of separation give differ- 
ent answers for the cooling flow fraction (for a review see 
ISunll2012|). which m ight be also affected by the sample bias 
l|Eckert et alll201ll ). 

In Table Table [7] and Table [7] we list all systems from 
both complete samples (B55 and HIFLUGCS) in order of 
right accession. Table [7] gives the X-ray peak and optical 
coordinates for the brightest cluster galaxies (BCGs), to- 
gether with the BCG names and the separation between the 
ICM's X-ray peak and the BCG center. The coordinates of 
the X-ray peak are well determined (within 1-2") for sys- 
tems that have a surface brightness peak; for the fainter 
and more diffuse systems the determination of the separa- 
tion becomes harder, and we mark these objects with an 
asterisk in Table [7] The optical coordinates for the BCGs 
are robustly determined (typically, within 1-2"), except in 
a few cases: A2163, A1689 (the optical image is of poor res- 
olution); A1367 which has 2 BCGs, NGC 3862 and NGC 
3842, both of them far from the X-ray center, but neither 
of them appear to be related to the central X-ray emission; 
and A2256 which has 2 comparable BCGs, NGC 6331 and 
UGC 10710, plus a few smaller ones (we choose NGC 6331). 
However, for some ACIS-I observations, the cluster cores are 
located in the CCD gaps, which affects the determination of 
the cluster centre (e.g., A399, A3158, A3562, A2255). 

2.2 Bubble Systems 

There are 22 and 26 bubble systems that are detected 
in the B55 and HIFLUGCS samples, respectively. There 
are, in total, 31 systems with bubbles. Among the 31 sys- 
tems, the bubble identification s for 24 systems are pub- 
lished ir^t!ieJiter^m^A2204 ((Sanders et al. 2009b), For- 



5), Fo: 

nax dShurkin et al.l I2008T) NGC 4636 (|Baldi et al.l 12009). 
NGC 5044 jDavid et"aT] l201ll ), and 20 systems in the 
iRaffertv et alj ( 20061 ) sample of clusters with bubbles. For 
NGC 507, III Zw 54, NGC 1550, A496, A3391, A1060 and 
A3532 we identified new bubbles by inspecting the Chan- 
dra data. For NGC 507, NGC 1550 and A496 more than 
one Chandra observation is available in the archive. In these 
cases, we used the ciao tool merge-all in order to add multi- 
ple observation for visual inspection of bubbles (see Table[T|. 

The bubbles were identified and measured using 
unsharp-mask images made from the exposure-corrected im- 
ages. For each system, a variety of different smoothing scales 
was investigated. The scales for which the contrast between 
the bubbles and the surroundings was the greatest were used 
for measurement. However, in some systems (e.g., Ill Zw 54, 
NGC 1550), the unsharp-mask image might reveal structure 
that is not related to bubbles. Thus, the unsharp-mask tool 
cannot be blindly used for bubble detection. Rather, deeper 
X-ray and radio images are needed to confirm the presence 
of marginally detected bubbles. As a result, some bubbles 
(e.g., Ill Zw 54, A1060, A3391, A3532) need to be further 
confirmed by deeper X-ray data. 

Since the bubbles are presumably created by radio 



lobes, we investigated the presence of a extended radio 
emission for these systems. The corona sy stems A3391 
and A3532 have large scale radio lobes l|Mauch et al.l 
120031; IVenturi et alJlioOlh much larger than the cavities we 
measured from the unsharp mask images. On the other 
hand, the corona syst em A1060 has very small radio lobes 
l|Lindblad et al.lll985l ) that fit well inside the cavities that 
we measured, and for the corona system III Zw 54 there 
is only a NVSS image available in which any lobe emission 
is not resolved. For the remaining syst ems, NGC 507 and 
NGC 1550 have l arge-scale radio lobes (jMurgia et al.ll201ll ; 
iDunn et al.ll201oT ). but A496 doesn't have any radio image 
in the literature besides the NVSS one. As in the case of 
corona systems A3391 and A3532, for NGC 507 the radio 
lobes are much larger than the cavity we measure, but for 
NGC 1550 the radio emission fills very well the cavity that 
we measure. 

As we underline above, for some systems there seems 
to be a discrepancy between the cavities sizes derived from 
unsharp-mask images and the extent of the radio lobe emis- 
sion. For the corona class systems, A3391 and A3532, it 
might be that the X-ray surface brightness contrast is too 
low and we are not able to measure the full extent of the 
cavities. In this case we might have underestimated the 
bubble power (enthalpy). For NGC 507, the unsharp mask 
shows a clear small depression in the west side of the clus- 
ter and this is the bubble that we list in Table [T] Addi- 
tionally, there app ear to be shallow deficits at the location 
of the radio lobes dParma et al. 19861; Ide Ruiter et al.lll986l ; 
iMurgia et al.ll201ll ; Giacintucci et al.ll2011bh . such as at the 
northern part of the western radio lobe, but deeper X-ray 
data is needed to confirm this. It might be that such cav- 
ities have a flattened "pancake" shape and are filled with 
entrained material, rather than the perfect void of material 
and spherical shape that we assumed when we did the calcu- 
lation, and as a result the contrast between the cavities and 
the ICM is weakene d. This is in agreement with results from 
IMurgia et al.l (120111 ) where it was found that this system has 
old lobes, and as a result they might have entrained mate- 
rial, expanded and leaked. As a result, if for NGC 507 one 
uses the radio lobes sizes and spherical geometry assump- 
tion, the bubble power might be overestimated by a large 
factor. 

Bubble radii and distances from the cluster centers are 
listed Table[T] Also, Table[T]lists the buoyancy time, the time 
required for the bubble to rise buoyantly from the cluster 
centre to its present location at its terminal velocit y, Vt (for 
a description of the analysis see lBirzan et al.ll2004r ): 



^buoy 

= R/v t ~ R^JSC/2gV, 



(1) 



where R is the projected radial distance from the cluster 
core to the bubble's center, S is the cross section of the 
bubble, V is the volume of the bubble and the grav itational 
acceleration g w 2<r 2 /R (jBinnev fc Tremaindll987l ). For the 
velocity dispersion ( a) the HyperLeda database] was used 
jPaturel et al.ll2003T) ; when no measurement of the veloc- 
ity d ispersion was available, a — 281 km s _1 (jBirzan et al.l 
|2004| ) was adopted, which was the average value. 
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Table 1. Bubble Properties 
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0.4 


0.000048 
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A3532* 


9.4 
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7.7 


7.7 


17.6 


0.36 


3.4 








S 


5.9 


5.9 


10.0 


0.16 


1.6 



"The sy stems for which we iden tified possible bubbles. The asterisk denotes the corona class 
systems iSun et alj|2007l ; ISunll2009h . 

6 For III Zw 54, A496, A3391, A1060, A3532, the time on source after re-processing the data (the 
same as in the Table[7]l. For NGC 507 and NGC 1550 the total time from multiple observations 
(Obs ID 5801 for NGC 1550). 

c For the velocity dispersion, the HyperLeda database was used ((Paturel et al.l[20 03l. When there 
were no measurements available, a = 281 km s" 1 was adopted IIBirzan et al.l l2004). 
^Projected semimajor axis of the bubble. 
e Projected semiminor axis of the bubble. 

^Projected radial distance from the cluster core to the bubble's center. 



2.3 Radio Properties 

Since the bubbles are thought to be the result of the inter- 
action between the central radio source and the ICM, it is 
of inter est to determine w hether a central radio source is 
present. iMittal et al.l l|2009l ) performed an radio analysis for 
the HIFLUGCS sample using literature information (e.g., 
NVSS, VLSS, FIRST etc) plus archived VLA radio data. As 
a result, for the HIFLUGCS s ample the presence of a cen- 
tral radio source is based on the lMittal et all |2009l ) analysis 
except for A754, A1060, A1650 and A1736. These sources 
show the presence of a central radio source coi ncident with 
the B CG location (see Table [7] for references'). IMittal et alj 
(2009) described a radio source as central if there is a radio 
source within 50h^ kpc of the X-ray peak, except A3562, 
A2142, A4038 and A3376 where a cut of I2h^ kpc was used. 
For the systems that are in B55 but not in the HIFLUGCS 
sample we established the presence of a central radio source 
based on information from the literature. We define a radio 
source as central if its location is coincident with the BCG 
location. Table [7] lists the presence of a radio source for our 
sample and the systems marked with an asterisk are the ones 
with central radio lobes. 

Also, Table [7J lists the 1.4 GHz monochromatic radio 
luminosities for the central radio sources, calculated as fol- 
lows: 



£l.4GHz = 47T_DlSi.4GHz(1 + z) 



(2) 



where a is the spectral index assumi ng S v ~ v a and 
Si.4GHz is the flux density at 1.4 GHz. IMittal et all (|2009l ) 
found that there is a good correlation between the total radio 
luminosity and the monochromatic radio luminosity (except 
for 2A 0335+096, A3376, MKW 3S and A4038). When avai l- 
able, the radio luminosities are from IMittal et al. | l|2009l ). 
Otherwise, they are computed using informatio n from the 
literature (e.g. FIRST, TGSS, NVSS surveys, lLal fc Raol 



2004 1; iBirzan et al.ll2004l; iLindblad et al.ll 19851; iGovoni et all 
20091 ; iBaum fc O'Deal Il99ll ; IVenturi et all 1200 ll ). In cases 



where a spectral index was not available, a value of -1 was 
adopted. For the systems with no detected central radio 
source, Table [7J lists the upper limits from the NVSS im- 
ages (the numbers without errors in the above table) or from 
IMittal et al.l (|2009T ) for the HIFLUGCS sample. 



2.4 Merging Activity 

Our systems range from giant ellipticals and poor clusters 
to massive clusters, and they often reside in active environ- 
ments with various degrees of merging activity. Since merg- 
ing activity is part of cluster formation and evolution and 
might be different for cooling flow versus non-cooling flow 
systems, we searched the literature for such information for 
all our systems. This information is listed in Table [7] How- 
ever, we do not make any distinction between major mergers 
and minor ones. 

Signatures of merging activity can be seen in optical, X- 
ray, and radio data. In optical data, merging may be inferred 
from the detection of substructures in the cluster, such as 
local peaks in the spatial distribution of galaxie s or localized 



veloc ity substr ucture. This is the ca se of A2142 (Owcr s et al. 

llLeonard et al. 



| 201l[) . A1367 l|Cortese et al.l 120041), A1689 
boill). A2063 and MKW3S llKrempec-Krygier fc Krygier 



20011 ). among many others 



Il999l ). A3395s (|Donnellv et al. 
(see Table 0. 

Furthermore, the X-ray data can reveal shocks that 
migh t be related to merger activity (e.g., A754. lHenrv et all 
2004), the presence of a si gnificant velocity gradient (e.g., 



A576, iDupke et all l2007al). detection of a low entropy 
gas (e.g., A3266. Finogucno v et al.l 120061), and detection 
of cold fronts |Owers et al. l2009al; iGhizzardi et ah 2010T) 



and sloshing activity (e.g., A1795, iMarkevitch fc Vikhlininl 



© 0000 RAS, MNRAS 000, 000-000 



The duty cycle of radio-mode feedback in complete samples of clusters 5 



l200lf). Cold fronts were first detected in Chandra data 
[ Markevitch et alj I2OO0I ; IVikhlinin et all 12001a! ') as contact 



discontinuities at the boundary of a cool, dense gas that 
is moving subsonically through a hotter, more diffuse sur- 
round ing gas (for a review, see Markevitch fc Vikhlininl 
120071 ). Using the B55 sample. iGhizzardi et al.l (|2010r ) found 

that a large fraction of systems ha ve cold fronts. 

Sloshing was introduced by iMarkevitch fc Vikhlininl 

l|200ll ) to explain cold fronts in the relaxed cluster 
A1795. Other exam ples of sloshing are found in A164 4 
l| Johnson et al.ll2010l ) and Ophiuchus dMillion et al.ll2010l ). 
among many others (|Markevitch fc Vikhlininl 120071) . Slosh- 
ing activity is thought to be a gravitational phenomenon, 
due to oscillations in the dark matter potential which 
migh t be set off by minor mergers ( Ascasibar & Markevitch 
20061 ). or hydrody namical in nat ure (jChurazov et al.ll2003l ; 



Fuiita et ail 12004 ). Furthermore. IZuHone et al.l (|201oh per 



formed simulations of sloshing activity, and they show that 
such activity might be responsible for quenching the cooling 
in some cooling flow clusters. 

Lastly, radio haloes and relics are indicators of merging 
activity. Hal oes and relics are thoug ht to be formed in ma- 
jor mergers l|van Weeren et al"1l201ll), and radio min i-haloes 
may be indicators of sloshing IjZuHone et alj|2011ah . In Ta- 
ble systems which have a radio mini-halo, radio halo or a 
radio relic are indicated. 



3 X-RAY DATA ANALYSIS 
3.1 General ICM Properties 

All systems were observed with the Chandra ACIS detec- 
tor in imaging mode, and the data were obtained from the 
Chandra Data Archive. Details of the observations are given 
in Table [3 

The Chandra data were reprocessed with CIAO 4.2 
using CALDB 4.2.1 and were corrected for known time- 
dependent gain and charge transfer inefficiency. Blank-sky 
background files, normalized to the count rate of the source 
image in the 10 — 12 keV band, were used for background 
subtraction^ 

The analysis of the X-ray data closely followed that of 
IRaffertv et all (HI). Briefly, X -ray spectra were extracted 
in circular annuli with ~ 3000 counts centered on the cen- 
troid of the cluster emission. For the fainter systems we used 
fewer counts (2000) for the inner region to allow derivation 
of the properties as close to the core as possible (e.g., Tri- 
angulum Aus, A2657, ZwCl 1215, III Zw 54, A2147, A1736, 
EXO 0423.4-0840, MKW 8). Also, in the case of the corona 
systems, due to the very small cores, we used for the in- 
ner region fewer than 3000 counts (2000 counts for A3558, 
1000 counts for A3376, and 500 counts for A3395s, 3C129.1, 
Coma). In some systems, due to the diffuse nature of the 
cluster emission or to substructure, the centroid was diffi- 
cult to identify precisely. These systems are noted in Ta- 
ble [7) Spectra and their associated weighted responses were 
made for each annulus using CIAO and were fit in XSPEC 
version 12.5.1 (|Arnaudlll996f ). 



See http: //asc .harvard. edu/ contrib/maxim/acisbg/ 



Gas temperatures and densities were found by depro- 
jecting the spectra with a single-temperature plasma model 
(MEKAL) with a foreground absorption model (WABS) 
using the PROJCT mixing model. We derived the cool- 
ing times using the deprojected densities and temperatures 
found above and the co oling curves from APEC spectral 
model lSmith et al.l (|200lh . The pressure in each annulus was 
calculated as p — nkT, where we have assumed an ideal gas 
and n m 2n a . Lastly, the entropy is defined as S — kTu^ 2 ^ 3 
and has units of keV cm 2 . 

To investigate whether or not the bubbles have enough 
energy to balance cooling, ones needs the luminosity of the 
cooling gas inside the cooling radius. Within the cooling 
radius, radiative energy losses must be replaced to prevent 
the cooling of large quantities of gas. We define the cooling 
radius as the radius within which the gas has a cooling time 
less than 7.7 x 10 9 yr, the look-back time to z = 1 for 
our adopted cosmology. To find the total luminosity inside 
the cooling radius, we performed the deprojection using a 
single-temperature model. Table[7]gives values of t coo i for all 
clusters, and r coo i and L x (< r coo i) for only our subsample 
of clusters that require heating (see Section r5.ll) , 

To derive densities as close to the core as possible, we 
additional ly use the onion peel ing deprojection method de- 
scribed in IRaffertv et al.l (|2008h inside the innermost annu- 
lus used for spectral deprojection, assuming that the tem- 
perature and abundance of the gas are constant in the in- 
ner region used in spectral deprojection and equal to the 
derived emission-weighted values. For a number of objects, 
constraints on the density profile inside the innermost an- 
nulus were not improved by this method (see Table [7]). For 
these cases, we report only the emission-weighted values 
from the innermost annulus. 



3.2 Cooling Time, Temperature Drop and the 
Thermal Stability Criterion 

Our goal is to understand the biases and selection effects 
in the detectability of current X-ray bubble samples and to 
place limits on the fraction of systems that can have enough 
bubble power to balance cooling. As a result, from these 
complete samples we need a subsample of systems in which 
feedback is expected to operate (i.e., cooling flow systems). 
The cooling time in the core is a basic selection criterion for 
cooling flows: if the cooling time of a system is smaller than 
its age, which is typically a significant fraction of the Hubble 
time, then one expects that this system needs heating to 
prevent cooling. 

In order to calculate the central cooling times we used 
archival Chandra data. The minimum radius at which we 
could derive reliable cooling times depends on the integra- 
tion time and proximity of the source (see Table [7| for the 
radius of the inner annulus used in the deprojection). Since 
we want to compare the cooling times for all the systems at 
a single physical radius, as close as possible to the nucleus, 
we computed the cooling time at 1 kpc using the surface 
brightness profiles to obtain densities (see Section l3.ll) . For 
some systems the extrapolation did not work well, as the 
surface brightness profile is too noisy or drops towards the 
centre or there is significant substructure that is inconsis- 
tent with the deprojection method. This was the case for: 
A401, A1367, A1736, A2256, A3667, A2319, A119, ZwCl 
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1215, A2657, A3562, A644 and A2255. For some other sys- 
tems, deprojection did not work because the surface bright- 
ness profile drops or flattens towards the centre or because 
they have X-ray-bright AGNs (such as Perseus) or large bub- 
bles. This is the case for A133, A478, A2597, Perseus, A2199, 
A4059 and Cygnus A. In these cases, the cooling time at 1 
kpc was extrapolated from the cooling time profile from the 
single-temperature model deprojection. This extrapolation 
was generally small, since for all these systems the inner bin 
has a radius below 10 kpc. For some systems there is a large 
difference between the cooling time from the innermost an- 
nulus (used in spectral deprojection) and the cooling time 
at a radius of 1 kpc. In some cases this is due to a very 
peaked surface brightness profile. This is the case for MKW 
8, A2634, A400, III Zw 54, A576, A1060, A3391, A2142, 
A3532, Coma, A3376, A3395s, and A3558. 

Additionally, we calculated the central temperature 
drop, since for a cooling flow cluster the temperature is ex- 
pected to drop towards the center. We calculated the drop 
as the ratio between the highest temperature in the profile 
and the temperature of the innermost annulus. Table [7] lists 
the temperature drop values for all systems. Some systems 
have no entry since the temperature profile increases towards 
the center, or the profile was too noisy and the temperature 
drop value was insignificant within errors. The calculated 
temperature drop depends on the size of the innermost an- 
nulus. However, the temperature typically varies slowly with 
radius, so it will not affect our estimates. Based on the tem- 
perature drop criterion only 2 systems with short cooling 
times at 1 kpc (less than 1 x 10 9 yr) do not show a temper- 
ature drop (A1651 and A2063). However, the temperature 
declines are probably related to the local virial temperature, 
so the cooler clusters, where the virial temperature of the 
central galaxy is closer to that of the cluster, have smaller 
central temperature declines. The temperature decline is not 
used in this paper as a criterion to separate the cooling flow 
systems from non-cooling flow systems, but just as an ad- 
ditional test after the systems where already separated by 
using a different criterion, such as cooling time. 

An alternative way to sel ect cooling flows is based on the 
therm al stability of the gas JVoit et al. 2008; Sh arma et al.l 
l2012f ). Using the sample of clusters. IVoit et al.l <|2008f > found 
that star formation and H-a emission (and hence cooling) 
seem to occur only if, at some location in the cluster, the 
following condition is met: 



??min = mm 



K{T)n e n H r< 



(3) 



where A(T) is t he cooling function calculated using APEC 
spectral model (|Smith et al.| [2001). and f c is the factor by 
which the magnetic field suppresses the conductivity below 
the Spitzer value. Assuming that the effective thermal con- 
ductivity can be expressed as a multiple, f c , of the Spitzer 
value, this parameter provides a measure of the stability of 
the gas to local cooling. For large values of Tjmm, thermal 
conduction overwhelms radiative cooling, preventing local 
cooling throughout the ICM. For small values of this param- 
eter, local cooling can run away, so that parts of the ICM 
may cool to low temperatures, resulting in the deposition of 
the cooled gas. If AGN are powered by the accretion of hot 
gas, then ?? m in should have little significance for the feedback 
process. However, if the AGN are fueled by cooled ICM and 



the effective conductivity is comparable to the Spitzer value, 
r7 m j n determine s the s ystems where fuel may be available. 

IVoit et all |2008t ) found that values of 7?min ^5 5 corre- 
spond approximately to an inner cooling time of 5 x 10 s yr. 
This agrees well with our findings (see Table [7] for the ?? m in 
values), the exeptions are A2065, A3391 and A400 which 
have ?7 m in < 5, but central cooling times of > 5 x 10 8 yr. 
However the l ast 2 systems ar e corona systems (with small 
cool cores, e.g. JSun et al .120071 ') with shallow Chandra obser- 
vations and we were not able to get closer to the mini-core. 



4 RESULTS FROM X-RAY AND OTHER 
OBSERVATIONS 

In this section, we examine how the central cooling time 
and rjmin relate to other properties of the cluster, such as 
the separation between the X-ray centroid and the optical 
centre (an indication of how relaxed a system is) and the 
central radio source luminosity. 



4.1 Cooling Time and Radio Luminosity 

Figure [1] left shows the central cooling time at 1 kpc versus 
the monochromatic 1.4 GHz luminosity of the central radio 
source. There is a separation between clusters, such that 
objects with radio luminosities above 2.5 x 10 30 erg s _1 
Hz -1 generally have cooling times at 1 kpc less than 5 x 10 8 
yr and have bubbles. There are no high radio luminosity 
ob jects with longer coolin g times. A similar plot is showed 



Cavagnolo et al.l |2008t ). but they were no able to probe 



the CF/NCF separation threshold since they used NVSS 
radio fluxes, and were not able to distingui sh between central 
radio source emission and halo emission. IsTml (2009) found 
a similar cutoff in radio luminosity for their strong cool-core 
syst ems (with cool i ng tim es below 1 Gyr). 

iRaffertv et al.l l|2008l ) studied star formation in a sample 
of 39 systems in different cooling stages and found a sharp 
threshold for the onset of star formation that occurs when 
the central cooling tim e falls below ~ 5 x 10 8 yr (see also 
ICavagnolo et al.ll2008T l. This number agrees well with our 
results shown in Figure[T]left. The only objects in our sample 
with a central cooling time > 5 x 10 8 yr and with a high radio 
luminosity (> 2.5x 10 30 erg s _1 Hz -1 ) are A3391, A400 and 
A1689. The first two are corona systems for which we were 
not able to probe the cooling time in the tiny cores. For 
A3391 we see bubbles in the existing Chandra image. The 
fact that there are no objects with high radio luminosity 
(> 2.5x 10 30 ergs -1 Hz -1 ) and long cooling times (> 5x10 s 
yr) indicates that the central AGN radio-mode heating in 
clusters is sensitive to the presence of gas with short cooling 
times. This finding fits well with the AGN feedback scenario. 

Also in Figure[T]we plotted in red the systems that have 
evidence in the literature for merger activity. We do not 
differentiate between a minor or major merger or sloshing 
(see Section ^. 4l for details). Figure [1] shows that mergers are 
present in both cooling-flow and non-cooling flow systems. 
We note, however, that the typical strength of the merger 
likely differs between the two classes. 
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4.2 Thermal Stability Parameter and Radio 
Luminosity 

IVoit et al.l l|2008l ) argued that a central cooling time of ~ 
5 x 10 8 yr is equivalent to a value of 5 for 7/ m i n , which agrees 
well with our results, shown in Figure [T] right, which plots 
rjmin versus the radio luminosity. A3391, A400 and A2065 
have central cooling times higher than ~ 5 x 10 8 yr and ?? m in 
smaller than 5. However the first two systems, as mentioned 
in the discussion of Figure [T] left, are corona systems with 
shallow Chandra observations and we were not able to get 
closer to the mini-core. On the other hand, A1689 has a 
central cooling time of ~ 5 x 10 8 yr, but ?7 m m much larger 
than 5. 

4.3 Cooling Time and X-ray Peak/Optical Core 
Separation 

Figure [2] shows a strong relationship between central cool- 
ing time and the separation between X-ray centroid and the 
optical centre (see Table 0). No object with central cool- 
ing time less than ~ 10 9 yr has a separation above 12 kpc, 
and almost no objects with a separation above 12 kpc (ex- 
ceptions A401, A119) have a cooling time below ~ 10 9 yr. 
For the systems with large separation between the optical 
centre and the X-ray centroid indicates that these systems 
are not dynamically settled that they are still going through 
significant merger activity. This trend is similar to the one 
in Figure \T\ left, but much tighter. The only system with a 
cooling time below ~ 10 yr that has a separation higher 
than 12 kpc is A4038. This system might have experienced 
major merger activity, which might have been cr eated the ra- 
dio r elic present on t he outskirts of t he clu ster (|Slee fc Rovl 
1 19981 ). Furthermore, iRaffertv etafl l|2008h shows that the 
systems with larger separation do not have star formation. 
This is another indicator that these systems do not need 
heating. 



5 THE AGN HEATING DUTY CYCLE 

5.1 Subsample of Clusters that Require Heating 

Our goal is to understand the biases and selection effects 
in the detection of X-ray bubbles, and to place limits on 
the fraction of systems that can have enough bubble power 
to balance cooling (these are the cooling-flow systems, CFs). 
However, in systems in which no cooling is expected to occur 
(i.e., non-cooling flows, NCFs), one also expects no heat- 
ing. Therefore, these systems should not be considered in 
our analysis. As a result, a subsample in which feedback is 
expected to operate must be defined from these complete 
samples. 

In Figure [1] (left panel) we found a clear separation be- 
tween CF and NCF systems. Furthermore, we found that 
the criteria using the thermal stability parameter and cool- 
ing time to separate CF and NCF systems give similar re- 
sults. Based on these results, we now define a subsample 
requiring heating based on the thermal stability parameter 
(»?min less than 5), since this criterium is more physically 
motivated. Additionally, based on the clear separation be- 
tween clusters seen in Figure [2] with a cut-off value of inner 
cooling time (at 1 kpc) of 1 x 10 9 yr, and a cut-off value of 




o.i L^j . . ,1 i I.. ,ii i J 

0.1 1.0 10.0 100.0 

Separation (kpc) 

Figure 2. Central cooling time (at 1 kpc) versus the separation 
between the X-ray peak and the BCG center. Filled symbols de- 
note the systems with bubbles and red colors denote the merging 
systems. Also plotted are the lines at 5 x 10 8 yr and 12 kpc. 

12 kpc for the projected separation between the X-ray peak 
and optical centroids, we further limited the subsample to 
systems that have such a separation of less than 12 kpc. 
In total, there are 49 systems (32 for B55 and 41 for HI- 
FLUGCS) that meet these criteria (rj m jn < 5 and separation 
< 12 kpc), 31 of which have detected bubbles (22 for B55, 26 
for HIFLUGCS). Table [2] lists the systems that meet these 
criteria and also lack detected bubbles. 

Finally, based on the central cooling time (at 1 kpc) 
and the separation (see Figure [2}, we also consider a some- 
what larger subsample: those systems with central cooling 
times below ~ 10 9 yr and separations smaller than 12 kpc. 
This subsample includes the previous main heating subsam- 
ple, plus 7 extra systems (see Table [2}. There are a total 
of 56 systems (39 for B55, 47 for HIFLUGCS) in this ex- 
tended subsample. Therefore, there are 16 systems in B55 
and 16 systems in HIFLUGCS (19 systems in total) that are 
excluded from the samples of systems that need heating. 

5.2 Simulations 

For the systems without bubbles from our extended subsam- 
ple (which have cooling times at 1 kpc smaller than ~ 10 9 
yr), we performed simulations in order to find the locations 
and the sizes of bubbles that can be present in the cluster 
and remain undetected. These systems are listed in Table [2] 

No attempt was made to simulate the effects of bright 
rims around the bubbles, since that would complicate the 
model greatly. When present, these can enhance bubble con- 
trast, making them more evident. As a result, if undetected 
bubbles had bright rims, we may be overestimating their 
possible energy content. Regardless of this, the modeling 
used here provides good limits on the energy that may be 
contained in undetected bubbles. 

We simulated a 3 dimensional cluster using a spher- 
ically symmetric j3 profile model for the emissivity (de- 
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Figure 1. Left: Cooling time (at 1 kpc) versus t he monochromatic 1,4 GHz radio luminosity for all 75 systems from both B55 and 
HIFLUGCS. The radio luminosity is taken from iMittal et al, (2009) were available, otherwise was computed using fluxes from the 
literature (see Table [7] for details and references). The filled symbols denote the 31 systems with bubbles (see Section 12,21 1, The red 
symbols denote the systems that appear to be undergoing some merger activity such as major mergers, minor mergers, or sloshing. This 
information is taken from the literature (see Table[7]for references). Over-plotted are the lines at 5 X 10 8 yr, 1 X 10 9 yr, 3 X 10 9 yr and 
2.5 X 10 30 erg s -1 Hz -1 (see text for details), right: Thermal stability parameter (see Section I3.2H versus the monochromatic 1.4 GHz 
radio luminosity. The symbols are the same as on the left panel. 



rived from the existing archive observations, see Table [2]). 
The parameters of the ft model were fixed to those derived 
from the azimuthally averaged surface brightness profiles. 
When a double-/? model was preferred over a single-/? model, 
the double-/? model was used in the simulation. A cube of 
500 X 500 x 500 elements was populated with the correspond- 
ing emissivities. 

For simplicity, spherical bubbles were simulated in this 
cube by setting the emissivities inside the bubbles to zero. 
However, it is important to notice that observationally it is 
unclear whether the X-ray cavities are devoid of X-ray emit- 
ting gas or instead are filled with hot, underdense gas. Then, 
by integrating the emissivity along lines of sight through the 
cube, we obtained the 2 dimensional surface brightness map. 
This map was used as input for the MARX simulatoi0 in 
order to obtain a simulated Chandra image. For MARX we 
used a spectrum extracted from the data, with the integra- 
tion time and position set as in the data. As in the case of 
the systems from Section [2~2l we used the ciao tool merge_all 
to add all the available observations for each system. This 
total integration time is used in the simulations (see Tabic 
[2] for details). 

We assume that the bubbles expand adiabatically into 
a /9-model atmosphere. Therefore, the radius of the bubble 
at at the distance r from the centre of the cluster is given 
by: 



Rb(r) =R h (0) 



1 + 



3/3/10 



(4) 



where Rh(0) is the bubble size at the cluster center, r c is the 
core radius and /? is the parameter of the beta model. To cal- 
culate the ages we assume that the bubbles rise buoyantly 
from the centre of the cluster to their current radius (giving 
ibuoy, see section T2.2I for details). In order to calculate the 
initial size, Rb(0), with which the bubbles are injected, we 
assumed: 4pV ~ Pcavt ~ Lxt, where p is the central pres- 
sure of the cluster, Lx is the bolometric X-ray luminosity 
inside r coo i (see Section [3]) and t is the time between the 
outbursts. We adopted a time between outbursts of 10 s yr, 
motivated by the observations of clu sters with multiple gen- 
erations of bubbles, such as Perseus (|Fabian et al . 2000). We 
ran multiple simulations for a cluster, with the locations of 
the bubbles and the angle in the plane of sky (<fi) random- 
ized, and with an angle from the line of sight (9) of 90° . The 
resulting images were inspected to determine if the bubbles 
remains visible for the whole of the 10 s yr outburst cycle. If 
a bubble become undetectable when its buoyancy age was 
less than 10 s yr, such a bubble may be present, but un- 
detectable, in the real cluster and the heating rate could 
be sufficient to balance cooling. In cases where the bubble 
would have remained visible for the whole cycle, the simu- 
lation was repeated using smaller values of P C av, until the 
bubble becomes undetectable in less than 10 8 yr. Under our 
assumptions, this value of P cav determines an upper limit on 
the possible AGN power in such systems. These systems are 
A3112, A1650, A1689, RXCJ 1504.1-0248, A2244, and Ophi- 
uchus. Therefore, these systems have upper limits less than 
the value required to balance cooling under our assumptions 
(see Table |2J. 



see http : / /www . space . mit . edu/ CXC/MARX/ 
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Table 2. Simulation Parameters 





t b 






ft 


J d 




f 


Kcsult 


System a 


(ks) 


(km s _1 ) 


(kpc) 








(kpc_) 






Main subsample: (M: 


inimum instability <5) and (Separation < 12 kpc) 






AWM 7 


180 


330 


6.8 


0.31 


121 


5.0 


4.7 


ix 


A400* 


21.5 


301 


135.6 


0.39 


5.4 


0.5 


4.3 


L x 


A3112 


109 




15.2 


0.43 


1091 


1.0 


10.5 


L x /5 


EXO 0423.4-0840 


25.9 




10.1 


0.44 


428 


0.17 


6.4 


Lx 


3C129.1* 


9.6 




66.3 


0.33 


11.5 


7.0 


6.4 


Lx 


A3376* 


60 


268 


136.1 


0.80 


5.7 


1.3 


5.7 


Lx 


A3395s* 


19.8 


335 


67.5 


0.38 


8.9 


1.4 


3.8 


Lx 


1V1JX VV 4 


30 


275 


2.5 


0.43 


455 


3.1 


A J 


-^X 


A1644 


70 




3.6 


0.31 


160 


1.0 


7.4 


Lx 


A1650 


228 




35.5 


0.42 


188 


0.41 


13.0 


Lx/5 


Coma (NGC4874)* 


9.5 


278 


140.7 


0.38 


35.8 


10.0 


6.0 


T 

^X 


A3558* 


11.2 


258 


37.5 


0.34 


72.3 


5.0 


11.0 


^X 


MKW 8* 


23.1 


217 


53.4 


0.31 


8.1 


0.5 


3.1 


1 


RXCJ 1504.1-0248 


53 




36.7 (145.8) 


0.75 (0.74) 


3319 (142) 


1.5 


18.4 


J /in 


A2065 


50 




53.0 


0.41 


65.4 


0.5 


10 


Lx 


Ophiuchus 


51 




4.3 


0.27 


416 


5.0 


12.0 


Lx/o 


A2589* 


93 


279 


3.0 (751.5) 


0.29 (21) 


108 (13.6) 


2.5 


8.5 


Lx 


A2634 (3C465)* 


49.5 


359 


0.05 


0.29 


397 


4.2 


4.1 


Lx 


Extra systems: (Cooling 


time at 1 kpc: 0.5-1 Gyr) and (Separation 


< 12 kpc) 




A576* 


25.6 




20.0 


0.30 


34.7 


1.0 


6.6 


Lx 


A1651 


8.9 




82.7 


0.51 


103 


0.29 


9.8 


Lx 


A1689 


194 




61.0 


0.54 


420 


0.54 


16.0 


Lx/5 


A3571 


31 


323 


10.9 (47.3) 


0.54 (0.33) 


77.7 (91) 


5.0 


5.6 


Lx 


A2063 


50 


222 


21.2 


0.35 


73 


1.6 


8.4 


Lx 


A2142 


44.6 




67.4 


0.48 


27 


1.0 


18.2 


Lx 


A2244 


56.7 




74.3 


0.55 


117 


2.3 


12.6 


Lx/5 



"Systems without bubbles in the existing Chandra images, but with a separation < 12 k pc and T? m ; n < 5 or a 
centr al cooling time (at lkpc) < 1 Gyr. The asterisks denote the corona class systems l lSun et alj|2007l : ISunl 
120091 ). 

''The total integration time from multiple observations when available: AWM 7 (Obs IDs 908, 11717, 12017, 
12018), A3112 (Obs IDs 2516, 2216, 6972, 7324, 7323), EXO 0423.4-0840 (Obs IDs 3970, 4183), A1644 (Obs 
IDs 2206, 7922), A1650 (Obs IDs 4178, 5822, 5823, 6356, 6357, 6358, 7242), A1689 (Obs IDs 540, 5004, 1663), 
RXCJ 1504.1-0248 (Obs IDs 4935, 5793), A2063 (Obs IDs 4187, 5795, 6262, 6263) and A2589 (Obs IDs 7190, 
3210, 6948, 7340). Otherwise, we listed the time on source after re-processing (the s ame as the ones in Table \7} . 
c Velocity dispersions were tak en from the HyperL eda database llPaturel et al.|[2003h . When none were available, 

a = 281 km s _1 was adopted ij Birzan e t alj |2004l) . 

d The normalisation of the isothermal surface brightness profile or /3-profile {Cavalicrc fc Fusco-Femianolll976t) 
in units of 10 — 9 ph s — 1 cm -2 pix -2 . 

e The background contribution for the isothermal surface brightness profile in units of 10 — 9 ph s — 1 cm -2 pix -2 . 
■^This is the bubble injection radius for the first run (L x , see Table[7]l; the injection radius for the actual run is 
5.8 kpc for RXCJ 1504.1-0248 (L x /10), 5.8 kpc for A1650 (L x /5), 4.7 for A3112 (L x /5), 7.1 for A1689 (L x /5), 
5.4 for Ophiuchus (L x /5), 5.64 kpc for A2244 (L x /2). 



5.3 Biases in Bubble Samples 

Samples of observed bubbles could suffer from biases due to 
detectability and projection effects. To assess the degree to 
which such biases might be present, we simulated an A4059- 
like cluster with two bubbles with a variety of sizes and 
positions. We then attempted to detect the bubbles by eye 
(without knowing their locations or sizes a priori). The mea- 
sured properties can then be compared to the intrinsic ones 
to understand possible biases. 

In total, 89 simulations were performed. Bubble radii 
were varied between 10 kpc and 60 kpc, and the distance 
from the core to the bubble was varied between 20 kpc 
and 80 kpc. Bubbles with radii greater than the bubble- 



to-core distance were not simulated. Figure [3] left shows a 
comparison between the measured bubble radii and the true 
ones. Approximately 44% of the simulated bubbles were un- 
detected. Of these, the majority had an axis between the 
bubble and core that was oriented near to the line of sight 
(6 < 45°) or were small bubbles at large distances from 
the core. Of those that were detected, the measured bubble 
radius is on average close to the true one. The greatest devi- 
ations from equality are generally for bubbles with smaller 
angles to the line of sight. 

A similar situation also applies to the comparison be- 
tween the measured distance to the core and the true one, 
shown in Figure [3] right. It is clear that most of the unde- 
tected bubbles lie at large distances from the core or, again, 
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''true (kpc) -Stmt (kpc) 

Figure 3. Left: Measured bubble radius versus the true radius. The color denotes the value of 6, the angle between the bubble-to-core 
axis and the line of sight. Undetected bubbles are shown at a measured radius of zero. Right: Measured distance from the core to the 
bubble versus the true distance. 




-Rbubble (kpc) 

Figure 4. Detection fraction as a function of the true bubble 
radius and the true (unprojected) distance from the bubble centre 
to the cluster core. The fraction has been summed over all angles 
between the bubble-to-core axis and the line of sight. 

have small values of 9. For detected bubbles, the average 
distance is close to the true one in most cases. These re- 
sults imply that when bubbles are detected, their properties 
are probably not grossly in error, particularly in an average 
sense. 

Lastly, we can estimate the detection rate for bubbles as 
a function of bubble radius and distance to the core by sum- 
ming over all angles of the bubble-to-core axis to the line of 
sight. This detection rate is shown in Figure[4] As expected, 
the detection fraction depends strongly on the bubble size 
and on the distance from the core. However, it is currently 
unclear as to which fraction of real bubbles are expected to 
fall into this parameter space. 

5.4 The Duty Cycle 

Figure [S] shows the bubble power (the heating rate) versus 
the bolometric X-ray luminosity (the cooling rate) within 
the cooling radius (at which t coo i < 7.7 x 10 9 yr) for our 



sample of cooling flows. The fraction of systems with de- 
tected bubbles for the main subsample is ~ 0.69 (22/32) for 
the B55 and w 0.63 (26/41) for HIFLUGCS. We note that 
this sample is based on two criteria: 77 m i n < 5 and X-ray-to- 
optical core separation < 12 kpc. If, however, the AGN are 
not fueled by cooled gas, then the r/ m i n criterion is not rele- 
vant to this analysis. For the larger samples based on cooling 
time cut-off the fraction is ~ 0.56 (22/39) for the B55 and 
« 0.55 (26/47) for HIFLUGCS. This implies a duty cycle in 
cooling flows of at least 55% to 69%. This duty cycle agrees 
well with previou s findings (|Dunn et al.ll2006l ; iNulsen et alj 
2009; iDong et al.ll2010T i and with the radio-loud fraction of 
BCGs in cooling flows (|Burnslll990h . It also agrees with the 
radio-loud f ractio n of ~ 30% found for massive galaxies by 
I Best et all (|2005l ). when adjusted for the fraction of such 
systems that are in cooling flows (~ 50%, see Section [5T]). 

From Figure \5\ we can conclude that most of the sys- 
tems that lack detected bubbles could have enough bubble 
power to balance cooling and remain undetected in exist- 
ing images. A few systems have upper limits on the bubble 
power below the 4pV line, implying that they cannot have 
sufficient bubble power to balance cooling under our assump- 
tions. However, one of our primary assumptions is that the 
bubbles are in the plane of the sky. As noted in Section [T] 
the orientation of the bubble axis relative to the line of sight 
has a strong effect on the detectability of the bubbles. If the 
angle from the bubble axis to the line of sight is < 30°, it 
is difficult to detect even large bubbles. As such an orien- 
tation is expected to occur in roughly 15% of all systems 
(assuming the jet axes are randomly oriented), it is possible 
that ~ 2-3 of the systems without detected bubbles in our 
sample have significant bubble power that is hidden. Un- 
fortunate orientation might therefore account for all of the 
systems for which the upper limit on the bubble power is in- 
sufficient to balance cooling. Therefore, with existing data, 
one cannot exclude the scenario that almost all cooling flow 
clusters have bubbles with enough power to balance cooling, 
implying a duty cycle for such activity of up to 100%. 

However, of the 7 intermediate cooling flow systems 
(with a central cooling time between ~ 5 - 10 x 10 8 yr), 
only A2063 shows extended radio emission. Thus, if the re- 
maining intermediate systems d o have bubbles, they m ust 
be "ghost", like those of A2597 (|McNamara et al.ll200lh . It 
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is more probable that these systems have previously had a 
cooling core that was heated by a merger, such as A2142, 
A3571, A1651, A2244, and corona system A576 (see Section 
I6.4|l . Therefore, these systems might not need bubbles in 
order to balance the internal heating. 

Also, it is important to mention that some of the sys- 
tems without bubbles are corona-class systems (denoted by 
the blue symbols in Figure |5J, which often have large radio 
lobes (all such systems in our sample except A576, MKW 
8, A2589 and A3558 have lobes) embedded in lower den- 
sity gas surrounding the bright X-ray core, which reduces 
the contrast of any bubbles and makes them harder to de- 
tect. Such systems may well have undetected bubbles. They 
include A400, A3376, A2634, A3395s, Coma, and 3C129.1. 

However, it is not fully understood what makes these 
systems different from larger cooling cores or whether the 
radio mode feedback cycle applies to them. Since the ra- 
dio lobes are inflated beyond the cooling region and the 
fraction of the outburst energy coupled to the corona is 
small, their feedback process may be very inefficient, and 
heating required to balance the internal cooling might b e 
supplied by SNe or by heat conduction (|Sun et al.ll2007h . 
Alternatively, radio heating through weak shocks may be 
important in these systems (for a discussion on the im- 
portan ce of weak shocks for AG N heating, see the re- 
view of lMcNamara fc Nulser]|2012l ). When the corona class 
systems are excluded, we find a duty cycle of 76% (= 
(22 - 3)/(32 - 7)) for the B55 main subsample and 77% 
(= (26 - 3)/(41 - 11)) for the HIFLUGCS main subsample. 

Lastly, we find that the fraction of systems with de- 
tected bubbles for both complete samples is very similar. 
The two subsamples overlap for 44 systems (see Section l2~Tj) . 
however the HIFLUGCS sample extends to lower-mass sys- 
tems. In general, our results are most relevant for clusters, 
since there are only 10 ellipticals and poor clusters in our 
sample (3 for B55 and 9 for HIFLUGCS), and 8 of them 
have bubbles. 



6 DISCUSSION 

6.1 Radio Luminosity Cut-Off 

Figure [T] left shows that there is a radio luminosity cut-off 
for massive clusters, such that systems with a central radio 
luminosity above 2.5 x 10 30 erg s _1 Hz -1 are found only in 
cooling flows (see Section [4.ip . This supports the idea that 
the radio-mode outbursts are self-regulating, with increased 
radio luminosity (and hence feedback power) seen in sources 
where cooling is likely occurring. 

We also note that the luminosity of central radio sources 
does not appear to be lower on average when there are no 
detected bubbles. In such systems, the radio source lacks 
bright lobes and is core dominated but is generally as bright 
as the systems with radio lobes and bubbles. This is the case 
of A1689, A1644, RXCJ 1504.1-0248, which represent ~ 10% 
of the high-luminosity objects (see Section ^. 3jl . A1689 is an 
interesting case, as it is not part of the main subsample since 
it has a thermal st ability parameter m uch higher than 5, and 
it has a radio halo (|Vacca et al.l201ll ) . This system may have 
been heated up by a major merger |Andersson fc Madeiskil 
l2004h that did not destroy the cool core. 



6.2 Duty Cycle of Radio Mode Feedback 

As discussed in Section f5. 41 our simulations imply an AGN 
heating duty cycle of at least ~ 60% and up to 100%. This 
result implies that cooling flow systems do not undergo long 
periods of time in which no significant cavities are present in 
their hot atmospheres. Therefore, it appears from our results 
that the atmospheres of cooling flows are undergoing almost 
non-stop energy injection through bubbles inflated by the 
central AGN. The exa ct nature of this heating rem ains un- 
clear (for a review, see lMcNamara fc Nulser]|2007f ). but our 
results clearly point to a connection between the presence 
of gas that is unstable to cooling (traced by low values of 
7?min) and the bubbles. This connection supports the idea 
that cold gas accretion, which might be induced by mergers 
or cooling of the ICM, fuels the outbursts that create the 
bubbles instead of hot accretion (i.e., Bondi accretion). 

The heating duty cycle that we have derived is distinct 
from the jet duty cycle, which is the fraction of time that the 
radio jet is active. Our results do not address the jet duty 
cycle, which, when traced only by the bubble population, 
could range from small values if the jets operate in short, 
intermittent bursts (which inflate the cavities and then stop) 
to values approaching unity if instead the jets are continuous 
(and the bubbles detach from them as they rise). However, 
we found that almost every source with cooling time below 
~ 5 x 10 8 yr has radio lobes. Our results should be useful as 
constraints on (or inputs to) simulations that include radio- 
mode AGN heating. 



6.3 Cooling Flow Systems without Bubbles 

Figure[T]left shows that the systems with bubbles have short 
cooling times ( < 5 x 10 s yr) and the majority of them 
have high radio luminosities (> 2.5x 10 30 erg s _1 Hz -1 ). 
However, as discussed in the previous section, some systems 
with short cooling t imes do not appear to have bubbles (e.g., 
iRaffertv et al.|[2008l ) . Here we further examine these systems 
to understand whether they are likely to have bubbles that 
are simply missed in observations. In particular, we examine 
their radio properties, star formation properties and merger 
activity. 

When the radio luminosity is considered, we find that 
many of the objects with high radio luminosity (> 2.5 x 
10 30 erg s _1 Hz -1 ) that do not show bubbles in the existin g 
Chandra images are "coronae" l|Sun et al.ll2007l ; |Sui3l2009h . 
In this category are A400, A3376, A2634, A3395s and Coma. 
The only coronae from this category for which we found 
bubbles are A3391 an d A3532. However, all these coronae 
have large radio lobes ([Owen ct al. 1985; Bagchi et al. 2006; 



Vcnturi et all 2001 



Hardcastle et alj|2005l : iBock et al.ll 19991 ; iBecker et aI1ll995l : 



The interaction between these radio 
lobes and the ICM is expected to create bubbles, as in A3391 
and A3532, but deeper X-ray images are needed to detect 
any bubbles embedded in the low density gas of these sys- 
tems. Another object with high radio luminosity, short cool- 
ing time and no visible bubbles in the existing Chandra im- 
age is EXO 0423.4-0840, which has a radio source at the 
centr e of the cluster with very small lobes (|Belsole et alj 
2005). Again, the interaction between these lobes and the 
ICM might create bubbles, but we are not able to see the 
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Figure 5. Bubble power versus the X-ray luminosity (inside the cooling region) for the subsamples of systems that require heating from 
B55 (left panel) and HIFLUGCS (right panel). The subsample of systems that require heating is listed in Table [2] (see Section 15.11 for 
details). Different symbols denote different subsamples: circles for the main subsample (r; m i n < 5 and separation < 12 kpc) and triangles 
for the extended sample (cooling time at 1 kpc between 0.5-1 Gyr and separation < 12 kpc). The filled symbols denote the systems 
with bubbles, and the blue symbols denote the corona-class systems, which are marked with an asterisk in Table [2] The diagonal lines 
indicate P C av = £x assuming pV, ApV or 16pV as the energy deposited. The limits on the 4pV line are not strictly upper limits, but 
simply mark the systems where the limit is at least this large. 



bubbles because the Chandra point spread function is too 
large to resolve them. 

The remaining category of objects that have high ra- 
dio luminosity, but lack bubbles, are massive clusters. This 
may be the case for A1 644, A3112 and RXC J 1504.1-0248. 
Except possibly A3112 dTakizawa et al.l2003h . these systems 
lack radio lobes. Therefore, these systems may well lack sub- 
stantial bubbles, with the possible exception of A3112, which 
has the shortest cooling time and some hints of radio lobes. 

These systems may be in a coo ling stage, between AGN 
outbursts. iBinnev fc Tabor] (|l995l ) postulated that there 
may be short bursts of nuclear activity alternating with peri- 
ods of quie scent cooling. Wh ile A3112 shows very little star 
formation (|fflcks et alj|2010h. RXCJ 1504. 1-0248 has a very 
high star formation rate (Ogrean et al. 2010; Raffcrty et al. 
2008) and a mini-halo l|Giacintucci et aHl2011aT ) most prob- 
ably related to the cooling flow activity, as in Perseus 
llGitti et al. 120021). A1644 also has high rates o f star for- 
mation (|Crawford et al.lll999l ; b'Dea et al.ll2010l ). Further- 
more, sloshing might pr ovide enough heat to balan ce the 
inner cooling of the gas (|ZuHone et al.ll2010L l2011bl ). This 
may be the case for A1644 which is thought to be under- 
going a merger between the main cluster and a subcluster 
i|Reiprich et aHl2004 Ijohnson et afl|201Ch ■ 

The cooling flow systems with a radio luminosity below 
2.5 x 10 30 erg s _1 Hz -1 are mostly ellipticals, groups or poor 
clusters (e.g., NGC 1550, NGC 4636, NGC 507, NGC 5044, 
Fornax, MKW 4, AWM 7, A262, A133) and most of them 
show bubbles in the existing Chandra images, except AWM 
7, wh ich is probably undergoing a merger ()Peletier et al.l 
Il99fj ). and MKW 4. Both AWM 7 and MKW 4 do not have 
any lobe radio emission, but they h ave star formation activ- 
ity (jMcNamara fc Q'Conneiilll989h . and thus they might be 
in a cooling stage, as proposed above for the more massive 
clusters like RXCJ 1504.1-0248. 

Another category of cooling flow objects with low radio 
luminosity are the corona class systems: A1060, III ZW 54, 
3C 129.1, A3558, A2589 and MKW 8. For A1060 and III 
Zw 54 we measured bubbles (see Section r2.2[) . and 3C 129.1 



has extended radio emission, and thus might have bubbles 
that are not visible in the existing shallow Chandra image. 
On the other hand, MKW 8, A2589 and A3558 do not have 
extended radio lobes (see Table [7}. 

There are also massive clusters with low radio luminos- 
ity: Ophiuchus, A1650 and A2065. The central galaxies in 
these clusters only show unresolved nuclear radio emission 
and no star formation, but A2065 and Ophiuchus might be 
in the midst of m ergers (see Table B . Also Ophiuchus has a 
radio mini- halo (|Govoni et all 12009] ; iMurgia et aHl2010ah. 
and shows evidence of sloshing in its core |N illion et al.l 
l2010h . which might be related (jZuHone etafl l2011al ). As 
in the case of RXCJ 1504.1-0248, the mini-halo in Ophi- 
uchus might be related to the previous cooling flow activity, 
and the s loshing has r e vived this old radio emission. Fur - 
thermore, lEckert et all {2008) and iNevalainen et all (|2009h 
show evidence of possible hard X-ray emission, which might 
be of non-thermal origin, caused by Compton scattering of 
cosmic microwave background radiation by the same popu- 
lation of relativistic electrons responsable for the mini-halo 
emission. As a result, Ophiuchus, A1650 and A2065 might 
be similar to the intermediate class systems, those with a 
cooling time between ~ 5 — 10 x 10 s yr, which presently are 
not expected to have bubbles since they are heated by major 
mergers, but may previously have had cooling flow activity 
and bubbles (see Section r6.4|l . 

Lastly, we note that small bubbles may be present in 
these systems but missed in current observations. Most of 
the B55 and HIFLUGCS systems are at moderate redshifts 
(around z — 0.05), with a few systems above a redshift 
of 0.1 (A1689 and RXC J 1504. 1-0248). Therefore, b ubbles 
such as those seen in M87 (e.g.. iFbrman et al.ll2005h would 
be too small to image. For example, if we extrapolate one 
of the inner cavities from M87, with a radius of 12 arcsec 
i|Forman et alj|2005h . to a higher redshift of 0.05 we will ob- 
tain a radius of 1 arcsec. As a result, this kind of bubble 
would be completely missed in current X-ray observations. 
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6.4 Cooling Flow Versus Non-Cooling Flow 
Systems 

The CF/NCF dichotomy has been studied a great deal re- 
cently using obse r vations (e.g.. [Sanderson et al.1120061. 
Chen et alTl2007l; ICavagnolo et al.l 120091 ; IPratt et all 



Hudson et alj|201oj) and simulations (e.g ., |Poole et al 



2009 



2009 



2008 



Burns et al.ll2008l : iMcCarthv et al.ll200Sl; iGuo fc Oh|[2009l) 



Using a statistically selected sample, Sanderson et al.l ( 20091 ) 
found that the dichotomy seems to be real and not an 
archival bias. There are two main underlying origins for 
the CF/NCF dichotomy. The first as sumes that the separa - 
tion occurs early d ue to pre-heatin g ( McCarthy ct al. 20081 ) 
through mergers (|Poole et all 120081; iBurns et al.l |2008|) o r 
from TeV gamma-rays from blazars (|Pfrommer et al.ll2012l ). 
The second interpretation assumes that the separation oc- 
curs late, such that a CF cluster can be destroyed due 
to merger (jRossetti fc Molendil |2010|; iRqssetti et alj l201ll ) 
or powerful AGN ou tburst (|Guo fc Oh! 120091 '). However. 
iPfrommer et al.l (|2012l ) argue that the AGN heating is likely 
insufficient to turn a CF into a NCF cluster, but the impact 
of AGN induced turbulence could resul t into a NCF atmo- 
sphere on time scale lar ger than 1 Gyr (|Parrish et al.ll2010l ; 
iRuszkowski fc Obll2010h . 

From Figures[T](left panel) and[2] we found that clusters 
are separated into two different groups based on the central 
cooling time, central source radio luminosity and the sep- 
aration between the X-ray peak and the optical centroid. 
Clusters with low central cooling times (< 5 x 10 8 yr) and 
smaller separation ( < 12 kpc) have a cooling flow activity 
and many of them show bubbles in their atmospheres. Clus- 
ters with longer cooling times (> 10 9 yr), low central radio 
luminosities (< 2.5x 10 30 erg s _1 Hz -1 ) and larger separa- 
tions (above 12 kpc) often show large haloes and relics likely 
due to merger shocks. 

For the cooling flow systems, most have visible bub- 
bles, or have extended radio emission from the central 
source (e.g., A3112), or appear to be in a cooling stage 
(e.g., A1644 and RXCJ 1504.1-0248). The systems with 
a cooling time between ~ 5 x 10 s yr and 10 9 yr (A576, 
A2244, A3571, A2142, A2065, A2063, A1651, A4038) do not 
have visible bubbles . This result agrees with the findings of 
Raffertv et all l|2008l ). since they are above the threshold for 
the onset of star formation (and hence significant cooling is 
probably not occurring). For the central radio source, none 
of them have more than point-like central radio emission 
(see Table 0, exc ept A2063 which ha s small radio lobes and 
possible cavities (|Kanov et al.ll200q ). Due to the fact that 
most of them are experiencing some sort of merger activity, 
diffuse emission might be also present in some of them, such 
as A2142, A3571, and A4038 (see Table[7]for references and 
M. Rossetti private communication for the detection of a 
radio halo in A2142). 

Among the 8 clusters with intermediate cooling times, 
only A2063 and A1651 lack the central temperature decrease 
typical of cooling flow clusters (see Tableland Section [3.2l) . 
Therefore, these systems share properties of both cooling 
flow and non cooling flow systems, and thus may represent 
an intermediate class of systems which are presently heated 
by a major merger but previously might have had a cooling 
flow activity and bubbles. However, one cann ot say whether 
they will evolve toward non-cooling flow as iRossetti et al.l 



|201ll) have suggested, or if they might restart the cycle of 
cooling and heating again at some point. 

Lastly, for the objects with central cooling times above 
~ 10 9 yr, none show any extended emission for the cen- 
tral radio source, except A2147 which shows small lobes 
(see Table 0, and they do not have visible bubbles A2 147 
might have been a corona s ystem which was destroyed by 
the merger dSun et al.ll2"007l ). However , this process ha s low 
probability even for corona systems ((Sun et al.|[2007h . All 
these systems with a central cooling time above ~ 10 9 yr 
might be experiencing majors mergers due to the continuing 
growth of large-scale structure. Therefore, they are classified 
as non-cooling flow systems for which one does not expect 
the presence of bubbles. 

We can conclude that cooling flow systems (including 
the intermediate systems) and non-cooling flow systems are 
well separated, expecially when we look at the central cool- 
ing time versus the separation between the X-ray peak and 
the optical core (see Figure [5J. Furthermore, the intermedi- 
ate systems separate from the cooling flow systems, with a 
central cooling time below ~ 5 x 10 s yr, if we look at the 
central radio luminosity (see Figure Q] left) and the presence 
of an extended central radio source (see Table [7}. The inter- 
mediate systems have a radio luminosity below 2.5 x 10 30 
erg s~ Hz -1 and they do not have extended central radio 
sources. However, from observations there is no conclusive 
way to know if a cluster like A2256 will become a cooling 
flow system (or if an object like A2065 will become a non 
cooling flow one). 

The merging NCF systems and the CF systems are 
quite separate in Figure [T] Thus any transitions between the 
two types must be rare and/or fast. However, a transition 
from CF to NCF caused by a major merger or a very power- 
ful AGN outburst was foun d to be improbable l|Poole et al.l 
2008]; IPfrommer et al.ll2012l ). Transitions in the opposite di- 
rection could be the result of gas cooling or the injection 
of low entropy gas during a merger. However, the speed 
of the former process is determined by the cooling time of 
the central gas, while it remains unclear whether the latter 
process is possible since the turbulence at the cluster cen- 
tre is sustained f or several gigayears after a cluster merger 
( Pau l et al . 2011). Further simulations are required to estab- 
lish which, if any, of these transitions are physically plausible 
and could be consistent with the data. 



7 CONCLUSIONS 

The duty cycle of radio mode feedback was determined for 
two overlapping, complete samples of cooling flow clusters: 
B55 and HIFLUGCS. All members of both samples have 
been observed with Chandra. As a result, these observations 
allow us with certainty to put limits on the non-detection 
of bubbles in some systems and to detect more bubbles in 
others. We found bubbles in 7 systems that were not pre- 
viously reported (NGC 507, III Zw 54, NGC 1550, A496, 
A3391, A1060 and A3532). 

We identified a total of 49 cooling flow systems (32 for 
B55 and 41 for HIFLUGCS) using two criteria: r/min < 5 
and X-ray to optical core separation < 12 kpc. However, 
we caution that r; m i n is relevant to AGN heating only if 
the AGN is fueled by cooled gas (as opposed to accretion 
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of hot gas, i.e., Bondi accretion). Our results imply that 
~ 60% of all clusters are cooling flows: ~ 57% for B55 and 
- 64% f or HIFLUGCS . Thes e numbers agree well with the 
result of iHudson etUI (|2010h for HIFLUGCS of 72%. If we 
exclude the corona systems (13 in total, 7 for B55, and 11 
for HIFLUGCS), the percentage of cooling flow systems is ~ 
45% (~ 45% for B55 and ~ 47% for HIFLUGCS). Out of 49 
cooling flow systems, 31 have bubbles (22 for B55, and 26 for 
HIFLUGCS). There are 7 intermediate cooling flow systems 
(7 for B55, and 6 for HIFLUGCS), and none of them have 
bubbles. Therefore, the duty cycle of radio-mode feedback 
is at least ~ 69% for B55 and ~ 63% for HIFLUGCS. If we 
exclude the corona class systems, we found a duty cycle of 
~ 76% for B55 and ~ 77% for HIFLUGCS. 

We found that the more luminous central radio sources 
(those with L1400 > 2.5 x 10 30 erg s _1 Hz -1 ) are only found 
in cooling flows (as are the X-ray bubbles). This result sup- 
ports the AGN feedback scenario, in which cooling from the 
hot atmosphere results in AGN outbursts (traced by their 
luminous radio emission) that in turn heat the gas and pre- 
vent periods of intense, long-duration cooling. 

However, some bubbles could be missed in existing ob- 
servations. For the systems without detected bubbles, we 
used simulations to limit the energy that could be present 
in undetected bubbles, finding that the power being injected 
by the central AGN into most of these systems could balance 
cooling, although the bubbles remain undetected. Therefore, 
with the existing data, one can not exclude the possibil- 
ity that all cooling flow clusters have bubbles with enough 
power to balance cooling, implying a duty cycle of up to 
~ 100%. This result needs to be investigated further, es- 
pecially using radio observations, to see if all cooling flow 
clusters have radio lobes that might create bubbles. 
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Table 3. X-ray and Optical Coordinates. 



X-Ray Core (J2000) a CDG Core (J2000) c Ar d 

System z RA DEC CDG Name 6 RA DEC (kpc) 



A85 


0.055 


00 41 50.2 


-09 18 11 


PGC 002501 


00 


41 


50.5 


-09 18 11 


4.8 


A119 


0.044 


00 56 15.9* 


-01 15 13* 


PGC 073498 


00 


56 


16.1 


-01 15 20 


6.7 


A133 


0.057 


01 02 41.7 


-21 52 48 


ESO 541-013 


01 


02 


41.7 


-21 52 55 


8.3 


NGC 507 


0.017 


01 23 39.9 


+33 15 21 


NGC 507* 


01 


23 


39.9 


+33 15 22 


0.2 


A262 


0.016 


01 52 46.3 


+36 09 12 


NGC 708 


01 


52 


46.5 


+36 09 07 


1.8 


AWM 7 


0.017 


02 54 27.6 


+41 34 43 


NGC 1129 


02 


54 


27.3 


+41 34 47 


1.4 


A400 


0.024 


02 57 41.7 


+06 01 36 


NGC 1128* 


02 


57 


41.7 


+06 01 35 


0.3 


A399 


0.072 


02 57 53.3* 


+ 13 01 31* 


UGC 02438 


02 


57 


53.2 


+ 13 01 52 


29.2 


A401 


0.074 


02 58 58.1* 


+ 13 35 00* 


UGC 02450 


02 


58 


57.9 


+ 13 34 59 


3.9 


A3112 


0.075 


03 17 57.6 


-44 14 17 


ESO 248-G 006 


03 


17 


57.7 


-44 14 18 


1.4 


Perseus 


0.018 


03 19 48.1 


+41 30 44 


NGC 1275 


03 


19 


48.2 


+41 30 42 


0.9 


NGC 1399 


0.005 


03 38 29.1 


-35 27 01 


NGC 1399* 


03 


38 


28.9 


-35 27 01 


0.2 


2A 0335+096 


0.035 


03 38 41.1 


+09 57 59 


PGC 013424 


03 


38 


40.6 


+09 58 12 


11.2 


III Zw 54 


0.029 


03 41 17.6 


+15 23 48 


2MASX J0341175+152348* 


03 


41 


17.5 


+ 15 23 48 
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04 
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04 
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04 
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04 
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06 
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0.7 


A3391 
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06 
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06 
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0.5 
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07 
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07 


47 


31.3 


-19 17 40 


2.1 


A644 
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09 


18 
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0.5 
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NGC 3311* 


10 


36 
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0.3 


A1367 


0.022 
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NGC 3862* 


11 


45 
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MKW 4 


0.020 
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NGC 4073 


12 


04 
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12 


17 
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M87 


0.004 
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12 


30 
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0.04 


NGC 4636 


0.003 
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+02 41 12 


NGC 4636 


12 


42 


49.8 
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0.3 


Centaurus 


0.011 
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NGC 4696 


12 


48 


49.3 
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1.0 
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2MASX J12571157-1724344 


12 
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3.3 
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2MASX J12572091-3021245* 


12 


57 
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12 
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12 
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13 
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13 


15 
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0.5 
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13 


26 
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A3558 
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13 


27 
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13 


33 
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14 


07 
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0.5 


MKW 8 
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14 


40 
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15 


04 
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3.5 
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0.077 
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+05 44 41 


PGC 054167 


15 


10 


56.1 


+05 44 41 


0.8 


A2052 


0.035 


15 16 44.4 


+07 01 17 


UGC 09799 


15 


16 


44.5 


+07 01 18 


0.9 


MKW 3s 


0.045 


15 21 51.7 


+07 42 23 


NGC 5920 


15 


21 


51.8 


+07 42 32 


7.8 


A2065 


0.073 


15 22 29.2 


+27 42 26 


PGC 054888 


15 


22 


28.9 


+27 42 45 


3.3 


A2063 


0.035 


15 23 05.3 


+08 36 37 


CGCG 077-097 


15 


23 


05.3 


+08 36 33 


3.2 


A2142 


0.091 


15 58 20.2 


+27 13 57 


2MASX J15582002+2714000 


15 


58 


20.0 


+27 14 00 


7.3 


A2147 


0.035 


16 02 15.4* 


+ 15 58 10* 


UGC 10143 


16 


02 


17.0 


+ 15 58 29 


20.6 


A2163 


0.203 


16 15 46.8* 


-06 08 39* 


2MASX J16154140-0609076* 


16 


15 


41.4 


-06 09 07 


283 


A2199 


0.030 


16 28 38.0 


+39 33 04 


NGC 6166 


16 


28 


38.5 


+39 33 06 


3.6 
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Table 3 — continued 



System 




z 


X-Ray Core 


(J2000) a 


CDG Name 6 


CDG Core (J2000) c 


Ar d 
(kpc) 


RA 


DEC 




RA 


DEC 


A2204 


o 
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IK 47 1 


4-0^ 34 34 


VT G>9 T1fi39 74-0^34 


16 


32 


47.0 


+05 34 33 


2 9 


r 1 Vi d vi (ti 1 1 1 1 m A 1 1 c 
_LI let lit U.1 Hill Jr\. Ho. 


o 


051 


it; 9Q 4* 


-64 21 11* 


F.^O 1 01 Ci 004 


16 


38 


18.1 
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43.2 


A2244 


o 
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1 7 09 49 fi 


4-34 03 38 

T^O 1 * UO OO 
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17 


02 
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9 2 


A99 c >fi 


o 
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17 8* 

J. I UO Jj.O 


4-78 3Q 00* 


MPP R331 * 


17 


1 


9fi 8 


+78 38 26 


117 


Onn 1 1 iflin t; 





.028 


17 12 27.7 


-23 22 07 


2MASX 117122774-2322108 


17 


12 


27.7 


-23 22 11 


2.2 


A2255 





.081 


17 12 41.3* 


+64 04 21* 


PGC 059830* 


17 


12 


34.9 


+64 04 15 


64.5 


A2319 





.056 


19 21 13.1* 


+43 56 11* 


CGCG 230-007 NED03 


19 


21 


10.0 


+43 56 44 


50.5 


Cygnus A 





.056 


19 59 28.3 


+40 44 02 


Cygnus A 


19 


59 


28.3 


+40 44 02 


0.5 


A3667 





.055 


20 12 31.6* 


-56 50 24* 


IC 4965 


20 


12 


27.3 


-56 49 36 


64.6 


Sersic 159/03 





.058 


23 13 58.7 


-42 43 31 


ESO 291-G 009 


23 


13 


58.6 


-42 43 39 


8.7 


A2589 





.041 


23 23 57.5 


+16 46 37 


NGC 7647 


23 


23 


57.4 


+16 46 38 


1.2 


A2597 





.085 


23 25 19.7 


-12 07 27 


PGC 071390 


23 


25 


19.8 


-12 07 26 


0.1 


A2634 





.031 


23 38 29.4 


+27 01 53 


NGC 7720* 


23 


38 


29.4 


+27 01 53 


0.1 


A2657 





.040 


23 44 55.9* 


+09 11 15* 


PGC 072804* 


23 


44 


57.4 


+09 11 34 


23.8 


A4038 





.030 


23 47 43.3 


-28 08 29 


IC 5358* 


23 


47 


45.0 


-45 08 27 


13.5 


A4059 





.048 


23 57 00.4 


-34 45 44 


ESO 349-G 010 


23 


57 


00.7 


-34 45 33 


10.5 



"The X-ray core positions from this work; the asterisks mark the uncertain core positions (due 



to the lack of the cooling core or to the presence of a confusing structure) . 

6 Notes on some systems (noted with an asterisk). NGC 507 — this system has a smaller satellite 
nearby. A400 — there are two cD galaxies in this system: A400-42, to the north, and A400-43 
to th e south. There is super-massive binary black hole at the centre of A400 jHudson et alj 
2006), we located the X-ray centre and the CDG centre on the northern AGN. NGC 1399 — 
this system has a satellite galaxy. Zw III 54 — a posible satellite galaxy is nearby. A3158 — this 
system has another cD galaxy in the centre part of the cluster, PGC 013652, and the location 
of the optical centre (PGC 013641) lies in the CCD gap. EXO 0422-086 — this system has a 
couple of satellite galaxies. A3266 and A3391 — the central cD galaxies in these systems appear 
to be interacting and have an extended corona. A3376 — this is not the central cD galaxy of 
the cluster, but the one that corresponds to the centroid of the X-ray emission. A576 — there 
is a another, similarly sized cD galaxy in the central part of the cluster, PGC020784. A1060 — 
another similarly sized cD galaxy is nearby, NGC 3309. A1367 — there is anoth er D galaxy, 
NGC 3842 located north o f the X-ray center, which has large scale radio lobes iBecker et alj 
ll995l ; IOwen fc Ledlowll 19971) , but neither N GC 3842 or NGC 3862 correspond to the X-ray peak. 
NGC 3862 hosts also a radio galaxy, 3C 264 llLara et al.lll999l ; ICondon et al.lll998ft . A3532 — this 
system has a satellite nearby. A1689 — there is another equally bright galaxy nearby; the optical 
image is of poor resolution. MKW 8 — there is another equally bright galaxy nearby. A2163 — 
the identification of the BCG is uncertain because the optical image is of poor resolution. A2256 
— the BCG identification is uncertain since there is another similarly sized cD galaxy in the field 
(UGC 10726), plus several smaller ones, such as CGCG 355-026 and MCG+13-12-014. Between 
the four cDs, only MCG+13-12-014 have an associated radio source (in NVSS). A2255 — there 
is another similarly sized cD galaxy nearby, PGC 059831. A2634 — this is a galaxy pair system; 
the pair is NGC 7720 NED01; A4038 — there is another smaller cD nearby, PGC 072437. 
c The CDG core positions are taken from the literature. 
d The separation between the optical centre and the X-ray center. 
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Table 4. X-ray Properties. 
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2 o+0.13 
u -^°-0.08 


3.4 


2+ ' 1 
u - z -o.i 


0.47 


4080l 4 ^ 


172 


A2029 


4977 


77.8 


o.Zo q 25 


n nc+0.01 

°- 06 -o.oi 


1.0 


06+ ' 01 
U ' UD -0.01 


3.84 


1049+5 


127 


A2052 


890 


36.8 




6+°' 5 

U.D_ 2 


2.8 


3+°' 2 


0.42 


9711 


80 


MKW 3s 


900 


55.3 


o 9 o+0.35 
°- zo -0.40 


6 +0 ' 4 

U.O_Q J 


3.3 


3+°' 2 


4.06 


104+J 


92 


A2065 


3182 


21.8 


1 89+ ' 54 
l.O»_ 44 


9 9+0.50 
z - z -0.5 


14.6 


7+°' 2 
u -'0.2 


2.68 


63+ 1 


63 


A2063 


6263 


16.8 




9 q+1.0 


7.2 


i-o±8:g 


18.9 


37+ 1 


68 


A2142 


5005 


44.6 


1 Q7+0.72 

1 - a '-i.o 


1 S+ ' 5 


6.7 


o.8l° : 4 


16.2 


7491^ 


126 


A2147 


3211 


17.4 


2 44+ 1 - 
-0.62 


10.0±\% 


20.0 


i-»±8:i 


6.0 
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Table 4 — continued 



System 


Obs. ID 


t a 


AT b 


'■cool 




icool (1 kpc) e 


77 ■ f 
'/mm 


ix 9 


r , h 

' cool 






(ks) 




(Gyr) 


(kpc) 


(Gyr) 




(10 42 erg s' 1 ) 


(kpc) 




IDOo 


C\7 A 

1 .4 




12 - 4 -5.5 


40. U 


p. fi+ 3 1 
6 - 6 -3.1 








A2199 


498 


15.9 


n ,^4-0 29 


r,o + 05 

0- 22 Iao5 


2.2 


. . -,+0 000 * 

0-117Io.oo5 


4.79 




106 


A2204 


499 


10.1 


ci +0.50 
,3 - oi -0.55 


n 1 k+0.03 
u - io -0.03 


1.6 


n nn+0.02 

0.09I 002 


0.44 


2024 2 g 


148 


Triangulum Aus. 


1227 


12.0 


n+1-2 

2-0-i.S 


a n+3.4 

6.9I 32 


12.5 


n + 1-2 
2 - 2 -l.l 


56.0 


•' " 




A2244 


4179 


56.7 


-, no + 0.17 
l - z6 -0.21 


-, n+0.3 
l -°-0.3 


5.7 


n K + 0.2 

- 5 -o.i 


12.1 


3I6I3 


113 


A2256 


2419 


9.2 




1S 7 +24.7 
!»•' _10.8 


27.4 






" ' , 




Ophiuchus 


3200 


50.5 


K ,r. + 0.53 

5 - 42 -0.37 


n oc+0.08 
°- 26 -0.09 


3.6 


n nnv+0.004 
°- 097 -0.005 


0.47 


324I 2 


83 


A2255 


894 


36.9 




15 -°-l 9 


33.0 




15.2 






A2319 


3231 


14.2 






25.5 










Cygnus A 


360 


33.7 


z - z '-0.33 


46+ ' 07 
u - 4D -0.06 


7.8 


n ir+o.08* 

u - iD -0.07 


2.1 


293+ 3 


92 


A3667 


889 


50.3 






24.0 










Sersic 159/03 


11758 


91.8 


2-30±g;f| 


n 99 +0.03 
u - zz -0.03 


2.5 


in+ 03 

u - lu -0.03 


0.64 


157+ 1 


111 


A2589 


7190 


53.4 




1 0+° 1 


1.8 


27+ - 04 
u - z '-0.04 


3.4 


2911 


62 


A2597 


7329 


53.3 


2.8+1? 


21 + 04 

u - zl -0.03 


3.2 


in+0.05* 
u ' iu -0.04 


0.64 


367+3 


119 


A2634 


4816 


48.5 


5 0+ 16 


n 77 +0.04 
u -' '-0.04 


8.36 


1O+0.09 
u - ±c> -0.08 


0.12 


6+ ' 4 
u -°-o .4 


17 


A2657 


4941 


15.7 




5.o±; : l 


12.7 










A4038 


4992 


29.9 


-, 91 +0.08 
1 - zl -0.07 


i-isl8:ii 


5.3 


1 0+ 1 ' 


9.7 






A4059 


5785 


92.1 


9 9n +0.08 
z - zu -0.08 




5.3 


n 0+0.1 * 
u -' 3 -0.1 


5.4 


85+1 


84 



"The time on source after reprocessing the data. 



6 The temperature drop, calculated as the ratio between the highest and the lowest temperature of the profile when the profile 
is rising upward smoothly. In some cases we neglected points that were outliers (e.g., in A576, A3376, A1689). For the systems 
which do not have a statistically significant temperature drop (e.g., A644) or their temperature profile is rising inward, there is 
no entry. 

c The cooling time of the innermost region (in Gyr). 

d The radius, in kpc, for the innermost region used for spectral deprojection. 

e The cooling time at 1 kpc derived using the deprojcctcd surface brightness profiles (in Gyr). In some cases the surface brightness 
deprojection did not work well (due, e.g., to the surface brightness profile dropping or flattening towards the center). In these 
cases, the cooling time at 1 kpc was extrapolated using the cooling time profile from the single temperature model deprojection. 
These systems are marked with an asterisk. For some systems a large extrapolation would be required. In these cases there is no 
entry for the cooling time at 1 kpc. 

•^Thermal stability parameter from lVoit et alj ll2008l) ; Inst. = min(kT/An e n[jr 2 ). For some of the systems the instability profile 
is still rising, therefore there is no minimum (e.g., A399, A3158, A754, A2163). See text for details. 

9 The X-ray luminosity within the cooling radius (determined from the sum of the deprojected fluxes for regions inside the cooling 
radius) for the subsample of systems that require heating in order to suppress the inner cooling (see Section l5,l| l, 
^Cooling radius, defined as the radius within which the gas has a cooling time less than 7.7 X 10 9 yr. 
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Table 5. Central Radio Sources, Bubbles, and Mergers 



System a 


Radio 6 


^1.4GHz C 




Bubbles 


Merger 


Sample^ 


AQt; 
rVoO 


yes^io,y±,yy^ 1 1 


no _i_ pi 1 7 w 1 pi 


-2 


ycsi^iuyj 


^r«c/^7zlS 79 190^ 

yes^o / j4o, i z, izy ) 


p. 1 ^ tttttt TTr'nt; 
jooo, xiir u uto 


Alio 

/\i iy 


no(z,4u,yo j 


o.y x iu 




no 


yes^OZj 


dck tTTTTT T TP PQ 




^roa{^ 1 iini 9R"\*t t 
yes^ii, i iu, izo ) 1 1 


1 s 4- n 1 v in - 
1.0 nz u.i x iu 


2 


yes^iuyj 


,.« e f1 1 n\ 
yes^iiuj 


nir J_j U utO 


TV no c;n7 

IN KjrVy ou / 


yes(ou,oi,yo,iuo ) 


(.10 it U.oO X IU 


-3 


ycs^ol,14oj 


ir( -io/'7/i i c\a ^ 
yes^ / 4, 1U4 j 


TTTTTT TTPPQ 
rllr 1j U \j O o 


A 9(^9 
iAZOZ 


yes i^i i ,zz j 


4.04 it U.U1 X IU 


-3 


Tree./' i ncA 
yes^iuy J 


i rod ( AQ\ 

yes (48 j 


bck TTTTTT TTPPQ 
JoOO, rllr Lt U b^S 


iAVV 1V1 ( 


nO(Z ) 


0.0 X IU 




no 


rvrnkoKln/ 1 07^ 

proDaDic^iu/ j 


rS00 


A Ann* 

/\4UU 


yes^oo,yo ^ 


1 A'\ 4- Pi 91 v 1 fl 
0.41 ^ U.Z1 X IU 


-1 


no 


yes (65) 


TTTTTT TTr i r ,( ^ 

nir 1j u vjOo 


A QOO 

Aoyy 


no(Z,y4 j 1 


l.o4 X IU 




no 


yes(y4, 110 J 


"Dec; TJTTTT TTi^r^C 

r>00, rllr LjULjL/o 


A401 


no(z j 1 


i O/i \/ in — 3 
l.y4 X IU 




no 


yes(y4,110J 


r>0o, rllrljUOL'O 


A ^1 1 9 


17»of Q1 1 98^* 

yesiyij izo ) 


1 79 4- Pt Pll 9 
1 . i Z it U.U1Z 




no 




Dec TTTTTT TTf^P^ 

r>oo, nir 1j u vjjV-/o 


Perseus 


t roc / 1 1 Q fi ~\ * 

ycs(ll,oo ) 


l./U ± U.U i 




ycs^iuy J 


ycs(Z4,40j 


rSOO 


WPP 1 QQQ 

in ±oyy 


yes(Z,iZ4 ( ) 


o qi j_ n ni v in 
z.yi it u.ui x iu 


-3 


yes (.1^4,) 




TTTTTT TTr^P^ 

nir 1j u vj o o 


z/\ uooo-tuyo 


„ oc /'l 1 1 1 Q 1 99 s !* 

yes^i i ,110, izz ^ 


1 f\« 4- n Pi 1 ^ v in 

l.UD It U.UO X IU 


-2 


ycs(^iuyj 


voc //io i /|9\ 

yes (40, i4z j 


R 1 ^ TTTTTT TTr^P^ 
r>00, nir LUuUo 


TTT 7™ KA* 
111 ZjW 04 


,, no p9 Q1 \ 

yes(z,yij 


q fit; 4- n in ^ in 
o.oO it u.iy x iu 


-3 


yes(l4d / ) 




TTTTTT TTPPQ 
nir 1j U \j L. O 


A ^1 ^,8 


yes^ooj 


q 1 1 j. n io v in 
y.ii ^ u.iz x iu 


-3 


no 


voo / 1 /i 1 \ 
yes(^i4i; 


RKtr TTTTTT TTr i r ,( ^ 

r>oo, nir luu o o 


A A 7S 


i roc / 1 1 ^ * 


qc; _|_ n 97 sy in 

o.yo ± Li.ii x iu 


-2 


yes^iuy j 




RKt; TTTTTT TTPPQ 

r>oo, nir Li u 


WPP i trca 


yes (04J 


^ 47 4- n v 1 n 

0.4( It U.OO X iu 


-4 


yes^oi , i4o j 


yes^i; 


TTTTTT TTf^PC- 

nir 1j u uLO 


TTVO 0A99 C\9.(\ 
L\iJ\\J U4ZZ-UOO 


yes (9)* 


a ^ 4-ni v in - 
4. i zn u. i x iu 


2 


no 




TTTFT TTf^P^ 

nir LUuL'D 


iAoZOO 


no(Z j 


i /ifi v i n — 3 
1.40 x iu 




no 


, rQC ,//|Q /1Q 191\ 

yes (40,40, Do, 1Z1 ) 


R^f; TTTTTT TTPPQ 
r>00, rilr LUVjL'O 




t roc f9^ 

yes(Zj 


9 n« j. n 11 v in 
z.yo zn u.ii x iu 


-2 


yes^i4dj 


^raafl^ A& 114 1 3^ 
yeS(00,4C5 3 1 14, lOO ) 


"R^ TTTTTT TTHP^ 
r>00, nlrLUvjL'O 


op i on i * 

o\_j izy.i 


yes( / 0, / 0, / y j 


1 nn _i_ n /i cc ^ in 
l.UU it U.40 X IU 


-4 


no 




r500 


Aoo / 


ycs(0,yi j ' 


1 m J_ n nfi w 1 n 
l.yi it U.UO X IU 


-1 


no 


ycs(0, 1U0 ) 


TTTTTT T TP PC 
nir Li U bljo 


A QQOI * 


ycs(o,00,yz,y / ) 


q nn J_ n /in 

v5.yy it u.4y 




yes(143J 


ycs(oo, loZ ) 


"Dec; TJTTTT TTi^f^C" 

r>00, nlrljULrL/O 


A ^Q^c* 

iT.ooyos 


yes(oj 


9 n^ 4- n no 
z.uo it u.uy 




no 


,, oc ('QQ 1 Q9\ 

yesi^oo, loz ^ 


rkk T4TTTT Tinn^i 
r>oo, nir ij u vjjV-/o 


AO/D 


yes (l,ylj 


0.4/ it U.Oo X IU 


-4 


no 


ycs(o0,4o, / o) 


per TJTTTT TTi^f^Q 
r>00, rilrLUbUo 


PT<f CI 074 ^ 1 Q1 


y cs V s .) 


^ PiQ 4- Pi 99 
O.OO It U.ZZ 




yes^uyj 




JjOO 


A £.A A 
iA.044 


n/~,/'9 s \ 

nO(Z J 


fi n ^ i r\— 4 

O.U X IU 




no 


! mo /l Q\ 

yes(10J 


rSOO 


A 7X.A 
r\ { 04 


, ^^i c./'7^^t 
y es v' u ; 


1 09 v in — 3 

l.UZ X IU 




no 


yes (4o, oz, 04, i u,oo j 


dc:c; TTTTTT TTr i r ,( ^ 

r>oo, nir l u uto 


Hydra A 


yes(ll,Z0,oU, loU J 


on i i n *"> 
oU.l it U.o 




yes^iuyj 




Dcr TTTTTT TTi^r^C 

rSOO, rilrLUbLo 


iAIUOU 


yes^oz; 


Q/19-tnQ/i \s in 

0.4Z it U.o4 X IU 


-5 


yes(14cij 


yes (61) 


Dec TTTTTT TTPPQ 
r>00, nlrLUbUo 


A 1 Qfi7 


no(z ) 


i a v i n — 4 

1.0 X IU 




no 


ycs(zo,4o J 


R^f; TTTTTT TTPPQ 
r>00, nlrLUbUo 


IVllv VV 4 


t roc /" Q 1 '\ 

yes^yi; 


91QJ_fl/ic; v 1 
Z.lo It U.40 X iu 


-4 


no 




TTTTTT TTr 1 ^^ 

nir l u uLo 


7„,pi 1 91 c; 
ZjWV_/1 1Z10 




A 7Q v 1 H — 3 

4. / y x iu 




no 




TTTTTT TTf^P^ 

nir l u \j o o 


Mo I 


ofi inn*!* 
ycs(ll,ZO,lUU ) 


7o _i_ n O/i w in 
0. 1 y it U.Z4 X IU 


-1 


ycs^iuy ) 


yes (llo) 


T2 r e 
rSOO 


WPP /IfiQfi 
IN^O 4000 


yes( / ,o4,oi ) 


9 S9 -U n m \/ in 

Z.oZ it U.UO X iu 


-4 


yes (7) 




TTTTTT TTPPQ 

nir j_j u \j ^ o 


Ccnta.ur ns 


yes^i i ,oy ) 


119 4 nq v 1 n 
1.1Z It U.UO X IU 


-1 


yesi.iuyj 


yes(45j 


dc:c; TTTTTT TTr^P^ 
r>oo, nir luvjL^o 


A 1 CKA A 


i roc / 9 

yes(zij 


a 7fi -U n 1 \s m 

4. (0 it U.10 X IU 


-2 


no 


(A Q (?7 1 1 1 \ 

ycs(4o,0 / ,111 j 


RKt; TTTTTT TTPPQ 

r>oo, nir ij u vjj^o 


A3532* 


yes(56,136)* t 


7.66 ± 0.30 x 10 


-1 


yes(143) 


ycs(136) 


B55 


A1650 


yes (55) 


7.7 ± 1.6 x 10" 


4 


no 




B55, HIFLUGCS 


A1651 


yes(2,91) 


1.32 ± 0.18 x 10 


-2 


no 




B55, HIFLUGCS 


Coma* 


ycs(l)t 


2.43 ± 0.01 x 10 


-2 


no 


yes(12,16,27,32,54) 


B55, HIFLUGCS 


A1689 


yes(l) 


7.77 ± 1.53 x 10 


-2 


no 


yes(5,81,134) 


B55 


NGC 5044 


yes(28,51)* 


5.95 ± 0.50 x 10 


-4 


ycs(28) 




HIFLUGCS 


A1736 


yes (4) 


2.92 ± 0.35 x 10 


-3 


no 




B55, HIFLUGCS 


A3558* 


yes(2,91)t 


2.11 ± 0.33 x 10 


-3 


no 


yes(48,135) 


B55, HIFLUGCS 


A3562 


yes(89)*t 


2.8 ± 0.5 x 10" 


4 


no 


yes(48,49,135) 


B55, HIFLUGCS 


A3571 


yes(137^ 


1.39 ± 0.06 x 10 


-3 


no 


yes(137) 


B55, HIFLUGCS 


A1795 


yes(ll, 15,46)* 


8.45 ± 0.24 x 10 


-1 


yes(109) 


yes(48,86,133) 


B55, HIFLUGCS 
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Table 5 - continued 

a 

oystcm 


IvaQlO 


c 

M.4GHz 




cubbies 


Merger 


d 

oamplc 


PKS 1404-267 


yes (66)* 


7.45 ± 0.26 x 10- 


2 


ycs(109) 




HIFLUGCS 


MKW 8* 


yes(l,91) 


4.02 ± 0.25 x 10" 


4 


no 




HIFLUGCS 


RXCJ 1504.1-0248 


yes(2,91)tt 


7.74 ± 0.27 x 10" 


1 


no 


ycs(50) 


HIFLUGCS 


A2029 


yes(l5,55,l3l)* t1 " 


1.06 ± 0.24 




yes(109) 


ycs(19) 


B55, HIFLUGCS 


A2052 


yes(ll,15,138)* 


1.56 ± 0.028 




yes(109) 




B55, HIFLUGCS 


MKW 3s 


yes(ll)* 


6.780 ± 0.001 x 10 


-2 


yes(109) 


yes(76) 


B55, HIFLUGCS 


A2065 


ycs(l) 


1.20 ± 0.22 x 10" 


2 


no 


yes(23,48) 


B55, HIFLUGCS 


A2063 


yes(91,99)* 


4.29 ± 0.28 x 10~ 


3 


no 


yes(76) 


B55, HIFLUGCS 


A2142 


yes(91)t 


4.60 ± 0.071 x 10~ 


-3 


no 


yes(59,84,102,103,133) 


B55, HIFLUGCS 


A2147 


yes(91,99)* 


4.37 ± 0.31 x 10" 


3 


no 


yes(77) 


B55, HIFLUGCS 


A2163 


no(2)t 


1.74 x 10" 2 




no 


yes(41,42,85,102) 


HIFLUGCS 


A2199 


yes(ll,47,53)*tt 


7.37 ± 0.16 x 10- 


1 


yes(109) 


yes(48,77) 


B55, HIFLUGCS 


A2204 


yes(119)* 


4.27 ± 0.15 x 10" 


1 


yes(119) 


ycs(120) 


B55, HIFLUGCS 


Triangulum Aus. 


yes (3) 


4.6 ± 0.8 x 10- 3 




no 




B55 


A2244 


yes(60,91) 


6.77 ± 0.82 x 10- 


3 


no 




B55, HIFLUGCS 


A2256 


no(2)t 


1.18 x 10- 3 




no 


yes(10,21,48,69,87,127) 


B55, HIFLUGCS 


Ophiuchus 


yes(2,4,93) tt 


5.17 ± 0.18 x 10- 


3 


no 


yes (90) 


B55 


A2255 


no(l)t 


2.34 x 10- 3 




no 


yes(88,54,108,117) 


B55, HIFLUGCS 


A2319 


no(2)t 


3.7 x 10- 4 




no 


ycs(48,54,96) 


B55 


Cygnus A 


yes(ll,17)* 


1092 ± 45 




yes(109) 




B55 


A3667 


no(2)t 


1.29 x 10- 3 




no 


yes(18,44,48,101, 102,115, 140) 


B55, HIFLUGCS 


Scrsic 159/03 


yes(ll)* 


1.92 ± 0.20 x 10- 


1 


ycs(109) 


probablc(29) 


HIFLUGCS 


A2589* 


no(2) 


5.83 x 10- 4 




no 




HIFLUGCS 


A2597 


ycs(ll,19,123)* 


3.32 ± 0.09 




ycs(109) 




B55, HIFLUGCS 


A2634* 


yes(15, 57,58)* 


1.730 ± 0.086 




no 




HIFLUGCS 


A2657 


no(2) 


5.49 x 10- 4 




no 




HIFLUGCS 


A4038* 


yes(125)t 


5.19 ± 0.08 x 10- 


3 


no 


probablc(14,125) 


B55,HIFLUGCS 
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Table 5 — continued 
System" Radio 6 



il.4GHz 



Bubbles 



Merger 



Sample 



A4059 



yes(ll,131)* 



6.82 ± 0.23 X 10" 



yes(109) yes(112) B55, HIFLUGCS 



References: (1) FIRST - Faint Im- 



ages of the Radio Sky at Twenty-cm l lBecker et alJll995T) ; (2) NVSS - The NRAO VLA Sky Survey at 1.4 GHz llCondon et alj 
Il998t) ; (3) SUMSS - The Sydney University Molonglo Sky Survey at 843 MHz llBock et al.l 119991); (4) TGSS (TIFR GMRT 
Sky Survey) by S . K. S irothi a, N. G. Ka n tharia , C. H. Ishwara-Chan d ra fc Gopal -Krishna; (5) lAndersson fc Madeiskil d2004l) 
(6) iBagchi et al.1 ll2006h: (7) iBaldi et , alj d2009ll; (81 iBaum & .O'Deal <ll99lh: (91 iBelsole et al.1 ll2005h: (10 1 iBerrington et al 
2002[) ; (lll| Birzan et al] 1120081 1; (121 iBravo-Alfaro et al.1 d2000h; (13) iBuote et all j2005l ); (14) iBurgett et alj d2004l); (15) iBurn; 



2005) : (21) IClarke fc Ensslinl (12006 
ICortese etaL ||2004|): (26) ICotton et al 
d2006h 
Dunn 



sjso) 

et al.1 



199C); (16) burns et alj dl994h: (17l]Carilli et al.l lll99lli': "(18) ICarretti et alj \201±; ( 1 9llClarke et alj fcool: (201 IClarke et al 



Clarke et all 120031): (23) IChatzikos et al.1 d2006l); (24) IChurazov et al l ||2003|); (25) 



de Ruiter et al 



(27) iDavis fc Mushotzkvl lll993ll; (281 iDavid et a .1 d201 ill; (291 



dl986h; (31) iDong etalJ d20ld); ( 32) 



(l2010l); (35) iDupke et al] d2007al); (36 ) 
et alj ||2005l); (40) IFeretti et al.1 Il999l 



(39) IFabian 

d2006T); (441 iFinoguenov et alj j201Ct>; (45) .. 
Ghizzardi et al.1 d2010l) ; (49) IGiacintucci et alj (|20Q. 
Giovannini et al. I dl993l); (53) iGiovannini et alj Jl99 



33) Donnelly et al 



de Plaa et alj 
ll200ll): (34 ) 



Donnelly et alj dl999T) . _ ... ._ 
Dupke et alj d2007bl); (3 7) iDurret etahl d2005l1; ( 38) IFabian et al.1 d200^1; 
(41) IFeretti et al.1 d200 lh ; (421 IFeretti et al.1 || 2004|1: (431 iFinoguenoy et alj 
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from lBirzan et al.l d2008h: for A1060 the 1.4 GHz flux was extrapolated from the 4.5 GHz one jLindblad et alJll985l) : for A1650 
the flux is from IGovoni et~al1 d2009h : for A1689 the flux is from the FIRST survey; fo r A1736 the TG S S sur vey at 150 MHz was 
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Galaxy Cluster Sample). See Section l2.1l for details. 
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